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Abstract 

Solvability of the rational quantum integrable systems related to exceptional 
root spaces G<i,Ft± is re-examined and for Eqj^ is established in the frame- 
work of a unified approach. It is shown the Hamiltonians take algebraic form 
being written in a certain Weyl-invariant variables. It is demonstrated that 
for each Hamiltonian the finite-dimensional invariant subspaces are made from 
polynomials and they form an infinite flag. A notion of minimal flag is intro- 
duced and minimal flag for each Hamiltonian is found. Corresponding eigen- 
values are calculated explicitly while the eigenfunctions can be computed by 
pure linear algebra means for arbitrary values of the coupling constants. The 
Hamiltonian of each model can be expressed in the algebraic form as a second 
degree polynomial in the generators of some infinite-dimensional but finitely- 
generated Lie algebra of differential operators, taken in a finite-dimensional 
representation. 
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1 Introduction 



Up to now the Hamiltonian Reduction Method which also is called the Projection 
Method P3 |2] provides a unique opportunity to construct non-trivial multidimen- 
sional completely integrable quantum (and classical) Hamiltonians. These Hamilto- 
nians are associated with root systems, they are related with the Laplace-Beltrami 
operators on symmetric spaces. Their rational (trigonometric) versions have re- 
markable properties: (i) their eigenvalues are known explicitly being at most the 
second degree polynomial in quantum numbers, (ii) any eigenfunction has a form 
of the ground state eigenfunction multiplied by a polynomial in the Cartesian (or 
exponential in Cartesian) coordinates. These two specific types of the Hamiltonians 
appear naturally in this approach - with rational and trigonometric potentials, cor- 
respondingly, characterized by the second-order poles in potentials on the boundaries 
of the configuration space. However, soon after a discovery of these Hamiltonians it 
became clear [2] there exists a straightforward generalization of these Hamiltonians 
to the case of arbitrary coupling constants without breaking any nice property. It 
led to a loss of immediate group theoretical interpretation. However, a property of 
solvability remained to hold. It gave a hint on existence of a more general formal- 
ism where above-mentioned Hamiltonians with arbitrary coupling constants appear 
naturally. An idea was to connect solvability with possible existence of an intrinsic 
hidden algebraic structure JUIIHE]. It turns out to be true. For arbitrary coupling 
constants these Hamiltonians admit an algebraic form (where polynomial coefficient 
functions occur in front of derivatives) and they are related with elements of the 
universal enveloping algebra of some algebras of differential operators acting in the 
space of invariants of the corresponding root space. Such an algebra was called the 
hidden algebra of the Hamiltonian. It was found that for all A n , B n , C n , D n , BC n 
rational and trigonometric models this algebra is the same (!) - it is the maximal 
affine subalgebra of the gZ n -algebra realized by the first order differential operators 
in R n and taken in symmetric representation jHE]. Thus, one can state that all 
these models are nothing but different appearances of a single model characterized 
by the hidden algebra gl n . Similar situation holds for the SUSY generalizations of 
above models - all of them turned out to be associated to the hidden superalgebra 
gl{n\n — 1), see [5]). 

However, one can naturally expect that the situation is drastically different for 
the Hamiltonians related to the exceptional algebras - each Hamiltonian is charac- 
terized by its own hidden algebra which is different for different Hamiltonians. A 
first indication stemmed from a study of the G2 rational and trigonometric models, 
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where the both models were characterized by the same hidden algebra, but this 
algebra turned out to be a certain infinite-dimensional, finitely-generated algebra of 
the differential operators g^ 2 \ which is a subalgebra of the algebra of the differen- 
tial operators on the real plane, C diff(2,R) [Hlllj. Later, it was shown that 
the similar situation holds for the rational and trigonometric F4 models jBJ: both 
models possess the same hidden algebra of the differential operators f^ 4 \ which is a 
subalgebra of the algebra of the differential operators on R 4 , C diff(4, R). 

For all previous studies a crucial point was to find a set of variables in which the 
Hamiltonian under investigation takes an algebraic form and reveal their meaning. 
A simple observation made in jl] was to consider the variables which incorporate 
all symmetries of the problem in hand. In particular, these symmetries contain (or 
sometimes even coincide with) the Weyl group of the associated root system. A 
natural idea is to take invariants of the fixed degrees of the Weyl group as variables. 
Many years ago V.I. Arnold [§] pointed out that the contravariant flat metric on the 
space of orbits of any Coxeter group written in terms of the polynomial invariants has 
polynomial matrix elements. It implies that the coefficient functions in front of the 
second derivatives in the Laplace-Beltrami operator are polynomials in invariants 
of the Weyl group. This result was rediscovered (and then generalized) later in 
111; 0, M an d |E] for A n , BC n , G 2 and F 4 algebras, respectively. It was shown 
that the algebraic structure persists for the whole Laplace-Beltrami operator: the 
coefficient functions in front of the first derivatives are polynomials as well. Further 
generalization was that similar statement is valid for the entire set of the rational 
Hamiltonians which are a combination of Laplace-Beltrami operator and a potential 
- after the gauge (similarity) transformation with the ground state eigenfunction 
the coefficient functions in front of the first derivatives remain polynomials 5 . In the 
present paper a certain universal prescription about a choice of coordinates is given. 

Recently, it was found that the existence of algebraic forms and knowledge of 
hidden algebraic structure of the above-mentioned Hamiltonians allows to consider 
their perturbations in a constructive way. One can develop an algebraic perturbation 
theory, where all corrections can be obtained by pure linear algebra means [T2*l 

cm. 

The goal of the present article is to carry out a study of the solvability of the 

5 Similar results we also obtained in above-mentioned works for the trigonometric A n , BC n , G2 
and -F4 Hamiltonians when the trigonometric Weyl invariants (which mean the Weyl invariants 
periodic in each variable) are used as coordinates. Even for elliptic A\ and BC n models the elliptic 
Weyl invariants allow to get the polynomials coefficients in front of derivatives (see ^U] an d d] 
for details) 
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rational models related to the exceptional root spaces. We introduce a general 
formalism which allows to study all these models on equal footing (Section 2). In 
Sections 3, 4 the solvability of the G2 and F4 models is re-examined in a new formal- 
ism, while in the Section 5-7 the solvability is established for £6,7,8 rational models. 
In Conclusion the results are summarized. Almost all results were obtained with the 
help of MAPLE 7 and MAPLE 8 programs together with the package COXETER 
created by J. Stembridge. 

2 Generalities 

Let us consider a quantum system described by rational Hamiltonian associated 
with a root system 1Z of algebra g of rank N: 



where a £ 1Z + are positive roots of the system 1Z which are vectors in R N , x = 
(xi, . . . ,xn) is a set Cartesian coordinates, \a\ 2 = Yli a \-> an d the scalar product 
(a ■ x) = ^ 1 ctkXk, 00 is a parameter. Coupling constants g\ a \ are assumed to be 
equal for roots of the same length. Hence for the A n case there is a single coupling 
constant, for the BC n case there are three coupling constants, etc. For some algebras 
(G2, E 6 , E 7 ) it is convenient to embed the roots into the vector space of higher 
dimension M 7V+n (n = 1 or 2). In general, the Hamiltonians of this type describe 
a quantum particle in multidimensional space, although for A n and G 2 cases these 
Hamiltonians allow another interpretation as the Hamiltonians describing many- 
body systems, while for B n , C n , D n they correspond to many-body systems on the 
space with mirror . 

Consider the spectral problem for the Hamiltonian 7i, 



where the configuration space is the Weyl chamber and \l/ should be from the cor- 
responding Hilbert space. Let us make a gauge rotation of the Hamiltonian taking 
the ground state eigenfunction \l>o as a gauge factor 




2 ' 



(2.1) 



(2.2) 



h = -2{^ Q {x))-\H-E^v{x) , 



(2.3) 
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where E is the ground state energy of the Hamiltonian (|2.1j) (see e.g. 2J), 

E = u(j+J2 9\a\j ■ (2.4) 

It should be mentioned that for the A n and G2 cases the energy E$ does include the 
ground state energy of center-of-mass motion. 

The ground state eigenfunction of ()2.1|) has a form 

tf = A ff exp(-|4 n) ) , (2.5) 

where 

and u\ a \ are defined through g\ a = v\ a \(v\ a \ — 1) and assumed to be equal for roots of 
the same length. For example, if all roots are of the same length as for A n case, we 
put v\ a \ = v. In the case of roots of two different lenghts as for G 2 case we denote 
v\ a \ = v for the short roots and V\ a \ = /i for the long roots. 4 is the invariant of 
the degree two (for definition see below). 

A new spectral problem arises 

hip(x) = —e<p(x) , (2.6) 

with a new spectral parameter e = 2(E — Eq). If in f!2.2|) the boundary condi- 
tion means normalizability of the eigenfunction \l/(x), then for (|2.6|) it requires the 
normalizability of <p(x) with the weight factor ^q(x). By construction, the lowest 
eigenvalue eo = and the lowest eigenfunction is <po = const. Our goal is to find, 
by a change of variables, an algebraic form of the operator h (if it exists). 

Definition: a linear differential operator with polynomial coefficients is called 
algebraic. 

In order to find these new variables we assume that they respect the symmetries 
of the Hamiltonian [4J. In our case it means the invariance under the group of 
automorphisms A g of the g root space. This group includes or coincides with the 
Weyl group W g . The algebraically independent invariant polynomials of the lowest 
possible degrees a generate the algebra S Aa of A 9 -invariant polynomials. The powers 
a are the degrees of the group W g . A particular form of these polynomials (denoted 



4 



below as t a ) can be found by averaging elementary monomials (u ■ x) a over some 
group orbit Q, 



(see e.g. ^H]), where x's are some formal variables. There exists a certain ambiguity 
in connecting these variables with the variables appearing in the Hamiltonian under 
study. It is worth mentioning that for any Lie-algebra g there exists a second degree 
invariant t\ (x) which does not depend on chosen orbit. It is namely this invariant 
which defines the exponent in the ground state eigenfunction (see above). One of 
natural connections is to identify x's with Cartesian coordinates. Later on we will 
use the expressions (12.7)1 as new variables in Hamiltonians under study and will call 
them the orbit variables. 

In general, making averaging over different orbits in (|2.7|) we obtain algebraically 
related invariants, except for a case when averaging over a specific orbit leads to 
vanishing results for certain invariant (s). Without loss of generality, these orbits 
can be discarded. In fact, the variables which were successfully used to solve the 
cases A n and BC n in [I] and [Sj, correspondingly, as well as for G 2 and F4 (see jHj 
and [Hj, correspondingly) can be easily obtained through the formulas ()2.7j) . 

The invariants of the fixed degrees are defined up to a polynomial in invari- 
ants of the lower degrees. This ambiguity plays an important role in obtaining the 
algebraic forms of Hamiltonians - these forms depend on the choice of variables, 
special combinations of variables only result in the simplest form which correspond 
to preservation of the minimal flag (see below). 

Now we introduce a notion of exact-solvability. Let us assume that the operator 
h possesses infinitely-many finite-dimensional invariant subspaces V n , n = 0, 1 . . ., 
which can be ordered forming an infinite flag 



(or, filtration) V. Therefore one can say that the operator h preserves the flag V. 
Definition [3] 

• An operator h which preserves an infinite flag of explicitly defined finite- 
dimensional spaces V is called exactly-solvable operator with flag V. We assume 
that the flag V is dense: between any two subsequent spaces there is no space 
of intermediate dimension which might belong to the flag. 




(2.7) 



V c Vi c V 2 c . . . c V n c . . . , 
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• If given h preserves several flags and among them there is a flag for which 
dimV n is maximal for any given n, such a flag is called minimal. 

Below we will deal with certain linear spaces of polynomials in several variables. 

Definition 

Consider the triangular linear space of polynomials in k variables 



where /'s are positive integer numbers and n is integer. Characteristic vector is a 
vector with components which are equal to the coefficients (weights) f\ in front of 



in k variables. The flag is defined using the induced grading. 

The characteristic vector is defined up to a multiplicative integer factor which 
we choose to be minimal. In most of the examples j\ = 1. Taking a sequence of 
the spaces characterized by growing integer numbers n we arrive at a flag which has 
<p(h,-,fk) ag g enera ti n g linear space. In the examples n takes consecutive integer 
values, n — 0, 1, 2, We call such a flag "pCft >•••>/&). 

All Hamiltonians we are going to study are of quite special type. All their flags 
of invariant subspaces, that we were able to find so far, are the flags of polynomials 
of the form ()2.8|) . Among these flags there always exists a minimal flag of a special 
form - comparing to other flags every fa takes its minimal value (!). For example, 
it was found that the minimal flag for G2 models (both rational and trigonometric) 
is V^ 1 ' 2 ' and the characteristic vector (|2.9|) is (1,2) 

Our final goal will be a search for minimal flags. It is worth to mention that 
one of the situations, when several flags of invariant subspaces can exist, occurs for 
the operator written in different variables while invariant subspaces remain to be 
polynomial one. Minimal flags have a remarkable property - they are preserved by 
corresponding trigonometric Hamiltonians if the latter are written in appropriate 
variables. 

A general strategy of our study is the following: (i) as a first step we gauge rotate 
the ground state eigenfunction, (ii) choosing a certain orbit we construct a particular 
set of variables which lead to an algebraic form of the gauge-rotated Hamiltonian, 



>p(f 



\,-,fk) 



<sf sf 2 . . .sf |0 < hVi + + • • • + fkPk < n) , 



(2.8) 
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(iii) exploiting ambiguity in definition of the invariants of the fixed degrees we search 
for variables preserving a minimal flag. 



3 The rational G2 model 

The rational G 2 model was introduced for the first time by Wolfes ^7j and later 
on was obtained in the Hamiltonian Reduction method 12 • This model allows 
an interpretation as a model of a three-identical particles with two- and three-body 
interactions. It was extensively studied in jS], [T2"j . 

The root system of the G 2 algebra is defined in 3-dimensional space with a 
constraint to the hyperplane y\ + y 2 + JJ3 = 0. Six positive roots are 

{e 2 - e x , e 3 - e 2 , -e a + e 3 , e x - 2e 2 + e 3 , -2e x + e 2 + e 3 , -e x - e 2 + 2e 3 } . 



The Hamiltonian of the rational G 2 model can be written in the form 
2 ^ L dx\ 



k=l 



k 



3 1 3 1 



where u; is a frequency parameter and g s = v{y — 1) > — |, ^ = — 1) > — | are 
coupling constants associated with the two-body and three-body interactions. The 
existence of two coupling constants reflects the fact that the root system contains 
two sets of roots, long and short: IZshort with roots of length 2 and lZi ong with roots of 
length 6. Parameter u in (jH.lj) is the only dimensional parameter in the Hamiltonian 
and the eigenvalues should be proportional to it. A connection to the root system 
with the model (jH.lj) appears when the center-of-mass coordinate is separated out 
and the relative motion is studied. Let us introduce the Perelomov coordinates of 
the relative motion 



yPerelomov = ^ _ ^ X , X = X X + X 2 + X 3 , (3.2) 

where X{ are the Cartesian coordinates of particles and X is the center-of-mass 
coordinate. Another convenient set of relative variables is given by the standard 
Jacobi coordinates 

Jacobi _ Jacobi _ Jacobi _ /o q\ 
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Both sets of the relative coordinates obey a condition y 1 + y 2 + 2/3 = 0. Thus 
the relative motion can be studied in three-dimensional y-sp&ce with the constraint 
2/i + 2/2 + 2/3 = 0. Transition from ()3.2j) to ()3.3|) corresponds to the interchange of 
two sets of roots - lZ s hort and 1Zi ong . One can identify the variables y either with 
yPereiomov QT yjacobi _ j^ ue ^ Q ^ e symmetry of the Hamiltonian with respect to the 
interchange of long and short roots a duality between the two sets of variables, 
yPereiomov anc j yjacobi^ a pp ears w j}[ demonstrated below. 

The ground state of relative motion is given by 

3 

<\V) = (AS r) (y)) i/ (A?°(y)) /1 exp{-|^^ } ) y l+ y 2 + y 3 = 0, 

k=l 

E = -u{l + 2v + 2/i) , (3.4) 
where A^\y) and A 2 \y) are Vandermonde determinants 

A^(x) = J] (a h ■ y) = f[(Vi ~ Vj) > 

T^-short J<* 

3 

A? (x) =l[( ak .y)= l[ ( Vl + y 3 -2y k ). 

Let us make a gauge rotation of the Hamiltonian (|3.1|) with the ground state 
eigenfunction (|3.4|) : 

fcg = -2(*«(y))- 1 (w2-^)*«(y), (3.5) 

and separate out the center-of-mass motion. 

Following the symmetries of the Hamiltonian (j3.1|) let us define the G 2 Weyl- 
invariant polynomials by averaging over the simplest orbit Q(e 2 — e\) , generated by 
(e 2 — ei). It has six elements: 

1 6 

k=l 

(cf. f)2.7|l ). where a = 2, 6 are the degrees of the G 2 -algebra and io k , k = 1, 2, ... 6 
are the orbit elements. Explicitly, 

4%) = 2(2/f + 2/1 + 2/12/2) , 
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4%) = 105^ 2 2 + 6&yly 2 + Myiyl + 22y« + 22y G 2 + lOOy^f + W&yfy* , 

Taking different orbits in the formula of averaging ()3.6|) we obtain different in- 
variants. They are related to each other as 

_ ,(n) 

fW = fW+^ n >(f< n )) s . (3.7) 

up to multiplicative factors. The general invariants of the lowest degree of the G 2 - 
algebra can be found through invariants obtained by averaging, they are algebraically 
related to tip\ In fact, an arbitrary set of invariants generating the S Ws > algebra can 
be found through a particular orbit invariants and has a form: 

t2 = tf\y) , 

h = tf\y) + A^\tf\y)f . (3.8) 

It can be shown that by taking t 2 ,6 as new variables in ()3.5|) we always arrive at 
the algebraic form of the gauge-rotated rational G 2 Hamiltonian for any value of the 
parameter One may assume that such a property should correspond to the flag 
with characteristic vector (1,3), which is invariant with respect to transformation 
(13. 8|) . The existence of this flag can be easily confirmed. 

A set of flags preserved by the Hamiltonian in coordinates ()3.8j) can be obtained 
by analyzing its action on monomials <f>p = t^t^ 2 labelled by vectors p = (p\,P2)- 
The monomial <j)f is mapped into a sum of monomials (p^j. ■ There are among di two 
zero vectors (0, 0) corresponding to terms —A&tid^ and —12u>t2dt 2 , which determine 
eigenvalues of the Hamiltonian. Other shift vectors are (—1, 0), (2, —1) and (5, —2). 
It means that the minimal flag for the general coordinates (J3.8|) is (2, 5). 

However, one can tune the value of the parameter to eliminate the term 
(5, —2) in the Hamiltonian. Remarkably there exist two such values of this param- 
eter, 

A {6) = -9/4 and A {6) = -11/4 , (3.9) 

resulting in smaller flag (1,2) (see below the equation (J3.13|) ). Following the defini- 
tion of Section 2 this flag p( 1,2 ) is minimal. 6 

6 It is remarkable that this flag is invariant under the action of the trigonometric Gi Hamiltonian, 
when this Hamiltonian is written in appropriate variables. 
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Making the substitution 



-4%: 

~ 4 



re = 4%)-?(4%)) 3 - (3-10) 



we arrive at the Hamiltonian 

h ( a? = 4r 2 ^ 2T2 + 24r 6 ^ 2T6 + 18t 2 2 t 6 «% T6 

- {Aujt 2 - 4[1 + 3(/i + i/)]} <9 T2 - [I2cur 6 - 9(1 + 2z/)r 2 2 ] d T6 . (3.11) 

This is an algebraic form of the rational G 2 model. It can be easily checked that the 
Hamiltonian (j3.11j) has infinitely-many finite-dimensional invariant subspaces 

Vt 2) = (rf rf | < Pl + 2p 2 < n) , n = 0, 1, . . . , (3.12) 

with the characteristic vector 

/ = (1,2) , (3.13) 

forming the minimal flag pt 1 ' 2 ) of the rational G 2 model which we denote p( G2 ) = 
(see a discussion above). 

The Hamiltonian (|3.11|) admits a representation in terms of a non-linear com- 
bination of the generators of some infinite-dimensional finitely-generated algebra 
g( 2 ) (f or definition and discussion see p|). It leads to the g^- Lie- algebraic form 
of the rational G 2 model. This form depends on a subset of the generators of the 
gl 2 x i? 3 -subalgebra of algebra g^ only and does not contain a raising generator 
of g( 2 \ The generators of g^ with excluded raising generator have infinitely-many 
common finite-dimensional spaces which coincide to the finite-dimensional invariant 
subspaces of the Hamiltonian ()3.11|) . 

We already studied the Hamiltonian ()3.5|) corresponding to the first choice of 
variables, when = —9/4 in ()3.8|) . Now we explore the second choice, when 
= -11/4. One can easily represent these variables in terms of the variables of 
the first choice 

f2{y)=tf\y)=r 2 , 

MV) = ~tf\y) + ^(4%)) 3 = -re + \rl , (3.14) 
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where y's are the Perelomov coordinates (|3.2j) . Making the change of variables (|3.14j) 
in (|3.11j) the Hamiltonian (J3.5j) takes the form 



4? = m A + 24f 6 ^ f6 + 18f 2 2 f 6 5| 6f6 

- {Aoof 2 - 4[1 + 3(/i + v)}} 3 f2 - [12WT6 - 9(1 + 2//)ff] 6> fe . (3.15) 

(cf. ()3.11Jl ). This is another algebraic form of the rational G 2 model, which is dual 
to ()3.11|) : the Hamiltonian ()3.15|) differs from f|3. 1 lj) by permutation v <-> /j, only. 
It reflects a symmetry between sets of short and long roots. It is evident that the 
Hamiltonian (|3.15|) preserves the flag V^ 1 ' 2 ' and admits a representation in terms of 
the generators of the algebra g( 2 > but now being written in the coordinates f's where 
v <-> ii. 

The duality between short and long roots also occurs when we make the identifi- 
cation of y variables either with the Perelomov coordinates ()3.2|) or with the Jacobi 
coordinates ()3.3|) . It appears in a relation between the r- variables (of the first choice 
of A^ 6 \ see (J3.9|) ) in the Perelomov coordinates and the f- variables (of the second 
choice of see (13 .9J1 ) in the Jacobi coordinates and visa versa, e.g. 

r 2 {y Jacobi ) = 2>f 2 {y Perelomov ) , T 6 {y Jacobi ) = 27f 6 {y Perelomov ) . (3.16) 

It is worth emphasizing that the coordinates f's f!3.14|) have a remarkable prop- 
erty: they are unique coordinates which can be 'trigonometrized' - they coincide 

. , , . . . ,. sin ax 

with the rational limit of certain trigonometric coordinates, hm a ^o — 

a 

in which the trigonometric G 2 model gets an algebraic form (see [6 ). The r 2 ^- 
coordinates do not have such a property 7 . 

It is interesting to analyze a transformation preserving the flag V^ G2 \ The most 
general polynomial transformation which preserve the linear space IS 

r 2 -> r 2 , 

r 6 -> r 6 + a {2 \\ , (3.17) 



7 It is worth noting that in the coordinates [ji , y/ro) the gauge-rotated rational G2 Hamiltonian 
(I3.5JI can be rewritten in terms of the generators of the maximal affine subalgebra of the gl($)- 
algebra [H] ■ Thus, this Hamiltonian possesses two different hidden algebras: the algebra acting 
on the configuration space parameterized by the r(f)-coordinates and the ^/(3)-algebra acting 
on the configuration space parameterized by the (r 2 , ^/tq) coordinates (see and Eqs.(2.6)-(2.7) 
therein) . 
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where is an arbitrary number of the dimension [co^ 1 ] for any n. It is evident 
that this transformation preserves the whole flag V^ G2 \ We find that there exist 
two algebraic operators which preserve the same flag V^ G2 \ However, the variables 
r's and f's in which these operators have being written are related one to each 
other through a transformation of the type (j3.17|) (see the relation (|3.14Jl ). Hence 
these two algebraic operators are non-equivalent. One suspects that the fact of the 
existence of two non-equivalent algebraic forms of the rational G 2 model (which look 
very much alike) reflects a certain intrinsic degeneracy of the model and should not 
hold for the trigonometric case. A study of the trigonometric G 2 model confirms 
that there exists a single g^-Lie-algebraic form. 

Thus, the operators (|3.11j) . (j3.15|) possess infinitely many finite-dimensional in- 
variant subspaces. These invariant subspaces coincide with the finite-dimensional 
representation spaces of the algebra g( 2 \ It is worth noting that for uj — v — fj, — 
both operators Hq coincide and represent the flat space Laplacian written in 
the (/^-Lie-algebraic form. It is evident that in this case there are no polynomial 
eigenfunctions in r(f)— coordinates. 

Both operators (ETTTj) . (l3~T5J) are triangular in the basis of monomials t^Tq 2 
(ff^fg 2 ). Therefore, the spectrum of (j3.11|) . (j3.15|) . h^ip = —2e(p, can be found 
explicitly and is equal to 



where n« are non- negative integers rix, n 2 = 0, 1, . . . (coefficients 2 and 6 are degrees 
of G 2 ). The spectrum does not depend on the coupling constants gi, g s (but the 
reference point for energy ()2.6|) does), it is equidistant and corresponds to the spec- 
trum of two harmonic oscillators with frequencies 2u and 600. Degeneracy of the 
spectrum is related to the number of partitions of integer number n to two weighted 
integers rii + 3n 2 . The spectrum of the original rational G 2 Hamiltonian (J3.1)) is 
E n = E + e n . It is worth noting that the Hamiltonian (j3.15J) (as well as (j3.11j0 pos- 
sesses a remarkable property: there exists a family of eigenfunctions which depend 
on the single variable T2(t2). These eigenfunctions are the associated Laguerre poly- 
nomials. This property allows to construct a quasi-exactly-solvable generalization 
of the rational G 2 model. It will be done elsewhere. 

It is worth mentioning that the boundaries of configuration space are determined 
by zeros of the ground state wave function (|3.4j) . In r- variables it is the algebraic 
curve 



e 



u(2ni + 6n 2 ) 





(3.19) 
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and in f-variables it has the same form 

■~\\2 



(AiA 2 (f ))' z = -27f 6 (f 6 - l/2ff) = . 



(3.20) 



In r 6 ,f 6 variables associated with the Perelomov ^-coordinates the product of dis- 
criminants looks especially simple 

(AiA 2 {y)f = 27r e (y)f 6 (y) , (3.21) 

where symmetry between short and long roots is seen explicitly. Also note ^Hj 



det 



dyk 



-i 2 



16 (AiA 2 ) 



(3.22) 



4 The rational Fa model 



The Hamiltonian of the rational F4 model written in the basis of the standard F± 
roots has the form (see |B] 8 ), 

w ^Ig ( _ c + ^ + s ,t(^ + _J_) (4 , } 



2 <— ' x 

i=i 



[ Xl + (-l)^X 2 + (-1)^X3 + (-1)^X4] 



2 ' 



where gfj = v(y — 1), g s — /i(/i — 1) are coupling constants related to sets of long 
and short roots, lZi ong and lZ s hort, correspondingly, and u is a frequency. Its ground 
state can be written as 



^ r) (x) = (A_A + r(A Arexp 



X i 



(4.2) 



i=i 



where 



A+A_ = JJ(a k -x) = Y[{ Xi + Xj ) Y[{xi - 



x 



3) ' 



n. 



long 



j<i 



j<i 



8 There is a misprint in [H] in the definition of coupling constants for the Hamiltonian (3.1): it 
should read gi = fi(fi - 1) > -1/4 , g = (l/2)u(v - 1) > -1/8 . 
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A A = Y[(a k x) 

T^short 

=rh n 

i=l v's=0,l 

The ground state energy is 

E = 2u{l + 6fi + 6u) . (4.4) 

Let us make gauge rotation of the Hamiltonian (|4.1|) with the ground state 
eigenfunction (|4.2j) as a gauge factor 

41 = -2(*W(x))- 1 (?4J - So)^(x) . (4.5) 

As new variables we take Weyl invariant polynomials of the lowest degrees of the 
group Wf a found by averaging elementary polynomials (u ■ x) a over the 24-element 
orbit Q generated by the root (ei + e<i) 

= ^ S> • *)° i a = 2, 6, 8, 12. (4.6) 

(cf. (12. 7j) ). In this case the powers a = 2,6,8,12 are the degrees of the group 
Wf 4 - As we mentioned in Section 2 the invariants of the fixed degrees a are defined 
ambiguously, up to some non-linear combinations of invariants of the lower degrees 

,(n) _^ ,(fi) 

4 n) -4 n) +A( 6 )(4 n) ) 3 , 

(Q) (Q) 4 (8)/.(0)x4 , ,(8),(n),(n) 
l 8 ~~ * r 8 + ^1 l r 2 J + ^2 r 2 r 6 > 

fg> - fg> + 4 12) (f< n) ) 6 + 4 12) (tf) 3 tf + Af } (4 n) ) 2 4 n) + Af\tf } ) 2 , (4.7) 

where J 4^ 6 ' 8,12 ' ) are arbitrary parameters. It can be shown that operator is alge- 
braic and it preserves a flag of polynomials for arbitrary values of these parameters. 
There are two flags (2, 6, 8, 12) and (2, 6, 8, 11) which are invariant under transfor- 
mations (|4.7p . They are similar to the flags (1, 3) and (2, 5) for the G2 case. 

Then we must search for a minimal flag. From a technical point of view it can 
be found by eliminating certain terms in the Hamiltonian, by choosing appropriate 



x x + (-l)-^x 2 + 



-V3 



Z3+(-l 



-1/4 



:r,i 



(4.3) 
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parameters A's. In jS] the minimal flag, denoted as V^ Fi \ is found. This flag is 
generated by the spaces of quasi-homogeneous polynomials 

7,(1,2,2,3) = (^VfVf rf | < Pl + 2p 3 + 2p A + 3j9 6 < n) , (4.8) 
with the characteristic vector / 

/= (1,2,2,3) , (4.9) 

which give rise to the flag 'p( 1 > 2 > 2 > 3 ). Hence, = T^ 1 ' 2 ' 2 ' 3 ' is the minimal flag for 

the rational _F 4 model. 

The characteristic vector ()4.9j) coincides with the highest root among short roots 
in the _F 4 root system written in the basis of simple roots 9 . Note that the most 
general polynomial transformation preserving the linear space *Pn ' ^ is 

Si — ► Si , 

s 2 — > s 2 + a 2 s\ + b 2 s 3 , 

S3 -> s 3 + a 3 s 2 + 6 3 s 2 , 

s 4 — ► s 4 + a 4 sf + b A sis 2 + c 4 sis 3 , (4-10) 



where a, 6, c are arbitrary numbers. 

Explicitly the set of variables preserving the minimal flag P^ 1 ' 2,2 ' 3 ) found in Ref. 
|H! is: 

r 2 = , 

12 6 12 v 2 ; 
8 80 8 30 2 6 48 V 2 ; 

r 12 = —4? - A(4«)) a 4n) + l(4 Q) ) 3 4 n) - — (4 n) ) 6 - — (4 n) ) 2 , 

720 12 288 v 2 ; 8 27 v 2 ; 6 1440 v 2 ' 1080 v 6 ; 

(4.11) 

In these variables the algebraic form of the gauge-rotated Hamiltonian ()4.5|) is: 

d 2 2, „ , d 2 , d 2 



4i "* — ^ r 2^2 + g(10 T 2 T 8 + T 2 2,i 6;^2 "T *V'2'12 ~T ^'8 '677^2 



9 We thank Victor Kac for this remark. The similar statement holds for A„-Calogero-Sutherland 
models where the flags are characterized by / = (1,1,.. .,1)0] and rational-trigonometric G2 
models where / = (1, 2) jjjj. However, it is not the case for the rational E$ model (see below). 
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d d d d d d 

+ 24r 6 — — + 32r 8 — — + 48r 12 - 



dr 2 <9r 6 dr 2 dr$ dr 2 dr 



12 



+ , r 2 2 r 8 + 6n 2 — — + 4 r 2 2 r 12 + 8r 8 2 — 7 — - 
3 dr 6 dr 8 dr 6 <9r 12 

d d d 2 
+ 4(3r 6 ri 2 + 2r 2 r 8 2 ) — h 6(2r 8 2 r 6 + r 2 r 8 Ti 2 ) 



<9r 8 9ri 2 <9rf 2 

- 4[ur 2 - 2(6z/ + 6/i + 1)]— - [12^r 6 - r 2 2 (4z/ + 2/i + 1)] — 

- 4[4^r 8 - r 6 (l + 3*/)] — - 4[6cur 12 - r 2 r 8 (2 + 3i/)]— (4.12) 

Cr 8 Cr 12 

The variables (|4.11|) are remarkable as they are the rational limits of certain trigono- 
metric variables in which the trigonometric _F 4 model takes an algebraic form [S|. 

However, the set of variables (|4.11j) does not exhaust all the possible sets of 
variables leading to the minimal flag ()4.9j) . Similarly to what happens with the 
rational G 2 model, there exists one more set of variables 

r 2 = t\ = r 2 , 

b 12 6 8 V 2 ; 24 2 

8 80 8 240 2 6 64 v 2 ; 8 4 2 6 192 2 ' 

~ __ 1 AO) 1 61 ^(Q)\a , 109 fr(»h2.(n) 847 f,(nh3 (n) , 109 / (nue 

ri2- 72Q ti2 + 17280^ 6 j + 5760 12 j 8 17280 1 2 j 6 + 3840 1 2 j 

13 1 1 

= -T12 + -To 2 r 8 + -r 2 r 6 T9 3 H r 9 6 , (4.13) 

8 2 8 6 32 6 2 2304 2 ' v ; 

(cf. (14.11)0 leading to an algebraic form of hy^ which preserves the minimal flag 
V {Fi) . The explicit expression for the Hamiltonian in the variables f|4.13|) gets the 
same form as ()4.12|) with v and /i exchanged: 

h%?{T)=h%?{T,v^v) (4.14) 

Clearly the Hamiltonian continues to be algebraic under the transformations (|4.1U|) . 

We were able to find the one-parametric algebra of differential operators for which 
Vi Fi) is the space of finite-dimensional irreducible representation (see Appendix B 
in Furthermore, the finite-dimensional representation spaces arise for different 



16 



integer values of the parameter. They form an infinite non-classical flag which 
coincides with V^ 4 ^ (|4.8|) . We call this algebra f^. Like the algebra introduced 
in jB] in relation to the G2 models, the algebra is infinite-dimensional yet finitely- 
generated. The rational F 4 Hamiltonian in either algebraic form (j4.12j) or (j4.14|) can 
be rewritten in terms of the generators of this algebra. 

The variables (j4.11|) and (j4.13|) can not be related by the transformation (j4.10j) . 
Therefore these two /(^-Lie-algebraic forms are non-equivalent from the point of 
view of the transformation (|4.1U|) . We interpret this fact by a certain intrinsic de- 
generacy of the rational _F 4 model. When the trigonometric F4 model is considered 
one can show that there exists the only one /^-Lie-algebraic form (up to a trans- 
formation (PTTUJO IB]. 

Thus, the operators ()4.12|) . ()4.14|) have infinitely many finite-dimensional in- 
variant subspaces. The finite-dimensional representations of the algebra f^. If 
uj = v = fi = 0, the operators coincide and become the flat space Lapla- 

cian written in the (/ 2 )-Lie-algebraic form, with no polynomial eigenfunctions in 
r(f)— coordinates. 

The operator (J4.12j) (as well as (|4.14j) ) is triangular in the basis of monomials 
rfV^rf 3 ^ ■ One can find the spectrum of (j4.12j) . h^ A <f = — 2ey2, explicitly 

eni,n 2 ,n 3 ,n4 = v(2ni + 6n 2 + 8n 3 + 12n 4 ) , (4.15) 

where = 0, 1, . . . are non-negative integers. Degeneracy of the spectrum is equal 
to the number of partitions of an integer number n to four weighted integers n\ + 
2>ri2 + 4n 3 + 6n^. The spectrum does not depend on the coupling constants gi, 
g s , it is equidistant and coincides (with different degeneracy) with the spectrum 
of the harmonic oscillator as well as with that of the rational -D4 model. Finally, 
the energies of the original rational _F 4 Hamiltonian ()4.1|) are E n = E + e n . It is 
worth noting that the Hamiltonian (J4.12|) possesses a remarkable property: there 
exists a family of eigenfunctions which depend on the single variable r 2 . These 
eigenfunctions are the associated Laguerre polynomials. As an illustration the first 
eigenfunctions are presented in the Appendix A. Again, this property allows to 
construct a quasi-exactly-solvable generalization of the rational F 4 model. It will be 
done elsewhere. 

Configuration space of the rational F 4 model (|4.1|) is defined by zeros of the 
ground state eigenf unction, i.e. by zeros of the pre-exponential factor in (|4.2j) . These 
zeros also define boundaries of the Weyl chamber (see e.g. (2])- The pre-exponential 
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factor (|4.2j) at v = /i = 2 can be written as a product of two factors. The first is 

(A + A_(r)) 2 = -192 7-^ + 256 rf , (4.16) 

it corresponds to the rational D4 model which occurs for g s = in (|4.1j) . The second 
one 

(A A(r)) 2 =-L(-192 rf 2 + 256 r| + 144 r 6 2 r 12 - 27 r 6 4 - 192 r 2 r 6 r 8 2 + 48 r 2 r 8 r 12 
4U96 

+ 30 r 2 2 r 2 r 8 - 12r 2 3 r 6 r 12 + ^r 3 r 3 + r 2 4 r 2 - \t\ r 6 r 8 + ^r 2 6 r 12 ) , 

(4.17) 

corresponds to the case of the degenerate F 4 model, gi = (which is equivalent to 
the -D4 model in dual variables, see [8 ). Thus, a boundary of the configuration space 
of the rational F4 model is the union of the algebraic hyper-surfaces (|4.16|) — (|4.17|) of 
degree 3 and 7 respectively, being the reduced algebraic hyper-surface of degree 10. 

In terms of the second set of variables one gets symmetric expressions 

( A + A_ (f )) 2 = -192 f 2 2 + 256 f 3 + 144 f 6 2 f 12 - 27 f 6 4 - 192 f 2 f 6 f 8 2 + 48 f 2 f 8 f 12 

+30 f 2 2 fl f 8 -12 r 2 3 f 6 fi 2 + ^ f 2 3 f 3 + f 2 4 f 8 2 - - f\ f 6 f 8 + - f 2 6 fi 2 , 

2 2 6 

(cf- dHISD ) and 

(AA (f)) 2 = ^(-192^ + 256^) , 

(cf.dHZD). 

In agreement with the general theory (see e.g. [T5] ) 

2 

= 4096(A+A_) 2 (A A) 2 . (4.18) 

To conclude a discussion of the rational _F 4 model, one can state that the ra- 
tional _F 4 model admits two non-equivalent algebraic and /( 4 )-Lie-algebraic forms, 
correspondingly. However, it does not manifest anything non-trivial. The existence 
of these two set of variables is related to a phenomenon of duality discussed in jH]: 
it can be easily checked that taking the variables 1)4.11)1 and substituting in ()4.6)1 
the variables rr's by z's, 

Zi j2 = xi±x 2 , 23,4 = x 3 ± x 4 (4.19) 
we get the variables ()4.13j) . 
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5 The rational Eq model 



The Hamiltonian of the rational Eq model is built using the root system of the 
Eq algebra. A convenient way to represent the Hamiltonian is to write it in an 
8— dimensional space {x\,X2, • • -x 8 } while imposing two constraints x 7 = xq, x 8 = 

^ 6 = -V) + ^ X> 2 + ^ 6 , (5.1) 



where 



5 

j<i=i 



+ 



1 



(/ > ■ . I / >"■ . ^ 2 ( / >-» . / y > . A 2 

1 



2 ' 



(5.2) 

5 

0, 1; z/j is even) 



is root-generated part of the potential with a coupling constant g = v[y — 1). The 
configuration space is given by the principal Eq Weyl chamber. 

In order to resolve the constraints, we introduce new variables 

Ui = 2-i j 1 — 1 • • • 5 

Uq = Xq + x 7 — x 8 , (using the constraint i/q = 3x 6 ) 

y-j = Xq — X7 , (using the constraint y 7 = 0) 

y 8 = Xq + x$ , (using the constraint y 8 = 0) 



(5.3) 



in which the Laplacian becomes 



d 2 



A (8) = A (5) + 3 + 2 

y dy£ 



Q2 Q2 Q2 



dy 2 dy£ dy 7 dy 8 



(5.4) 



while the potential part of (|5.1j) depends on {yi . . . y e } only: 



V 



LO' 

~2 



i=l 
5 



vl 



j<i=l 



4K^-EL(-iNi 



2 • 



(5.5) 
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In this formalism imposing constraints implies a study of eigenfunctions having no 
dependence on y 7 , y 8 . Hence, y 7i8 -dependent part of the Laplacian standing in square 
brackets in (|5.4|) can be dropped off. 

The ground state eigenfunction has a form 



,(5) A (5)VA« „ ' !r '- 



■$ Q = (AX > A™) v A v E6 e 2 I > , E = 3w(l + 12z/) (5.6) 

where 

aL 5) = n(w±yi) 

J<i=l 




with g = v[y — 1). 

In order to find variables leading to the algebraic form of gauge-rotated Hamil- 
tonian, 

h ( £(yi...y 6 ) = -2S> - 1 (H E6 -E )(y 1 ...y 6 )S> , (5.7) 

let us define a basis in the form of the Weyl-invariant polynomials averaged over the 
27-element orbit generated by the vector e^, 

27 

t<?> = J^i^ • y)\ oo k G fi(ea) , (5.8) 



k=l 



(cf. (|2.7jl ). where a = 2,5,6,8,9,12 are the degrees of the E 6 Weyl group and 
LUk, k = 1, 2, ... 27 are the orbit elements. The orbit variables y in (J5.8|) are identified 
with variables y\ . . . y§ in (J5.3j) and in (|5.8|) . The Weyl-invariant polynomials of the 
fixed degree are defined ambiguously 



.(H) _ (O) 
l, 2 fc 2 
t (0) _^ t (ct) 



b 2 

An) An) As) An) An) ,(8),,(n)s 4 

r 8 ~~ *■ r 8 + A l l 2 l 6 ' A 2 \ l 2 1 J 

_> t («) + ^)(tf) ) 2 tf , (5.9) 
,(n) An) 4 (i2).(n)/.(n)N2 , 4 (i2)/.(n)x 2 .(n) , ,(12) ,,(fi)s 3 ,(o) 

r 12 — *■ r 12 + A l r 2 l r 5 J + ^2 l r 2 ) l 8 + ^3 l r 2 J r 6 
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+4 12) (4 n) ) 6 +4 12) (4 n) ) 2 - 

where Af' 8 ' 9 ' 12 ^ are parameters. For general values of these parameters one gets 
the algebraic Hamiltonian preserving the flag 2, 5, 6, 8, 9, 12 and even smaller one 
1,2,3,4,4,6. Our goal, as before, is to tune these parameters in such a way that 
the Hamiltonian in its algebraic form preserves the minimal flag. Quite cumbersome 
analysis results to a one-parametric set of variables 

n - —t {n) 

r 6 = 3456 4 n) -2A(tf ] f , 

r 8 = 2 ^tf ) +4&{tfy- h ^tfhf\ (5.10) 



Tg 



663552 4 Q) 27648 4^4^ 2 



7 

^ = — Jl ^- t w _ rrp + 294912 ti l> (4 >,J ) 8 + A?*>t$" (CO 



95551488 . (n) 5568 (n) 6 0Q , Q10 .(n) ,.(n)x 3 , 4 (i2).(n) / .(n)x 2 



-1658880 (4 C) ) 2 4 n) -^^(4 n) ) 2 , 

where is a parameter, leading to the minimal flag, which we denote V^ E(> \ This 
flag is spanned by the spaces 

^(1,1,2,2,2,3) = ( r |vrTf 6 rr-r| >9 -rf 2 12 ] <p 2 +p 5 +2p 6 +2p 8 +2p 9 +3p 12 < n> , (5.11) 

with the characteristic vector, 

/ = (1,1,2,2,2,3) . (5.12) 

(cf. fl3.13j) and (|4.9j) ). which coincides with the Eq highest root, confirming Kac's 
conjecture. 

The most general polynomial transformation which preserves a linear space 
pC^e) _ -pO-, i,2,2, 2,3) ^ an _y n . g q ^ ^ e f orm 

«5 -> s 5 + a 5 s 2 , 
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s 6 — > s 6 + «6,l s 2 + a 6,2<S 2 + &6,l s 8 + be t2 Sg + ^6,3 s 2-55 , 

s 8 -> s 8 + a 8)1 S2 + a 8j2 S5 + 6 8i is 6 + 6 8i2 s 9 + b 8t3 s 2 s 5 , 

s 9 ~~ *• s 9 + Ct8,l s 2 + a 8,2-55 + &8,l s 6 + ^8,2-Ss + &8,3 S 2S5 , 

3 3 l 2 l 2 

s 12 ~~ * s 12 + &12,l s 2 + a 12,2'S5 + 012,1 S 2 S 5 + Ol2,2S2S 5 + C^l^^ + C12, 2S 2 S 8 

+Cl2,3S2S 9 + C?12,lS5S6 + dl2,2«5S8 + ^12,3^559 , (5.13) 

where {a, 6, c, c?} are arbitrary numbers. Surprisingly, there is an overlap of non- 

_____ _ _ M ty\ 

linear transformations (|5.1()j) and ()5.13j) . Namely, a variation of the parameter A\ 
in corresponds to varying the parameter 612,2 in the transformation ()5.13j) . 

Therefore there exists a non-trivial one-parametric set of invariants of fixed degrees 
leading to an algebraic form of the Hamiltonian h^ e and simultaneously preserving 

a minimal flag ! Thus, the parameter A\ can be chosen following our convenience. 
(12) 

We set A\ =0, which makes the coefficient functions in the algebraic form of the 
Hamiltonian (see below the polynomials of lowest degree. However, it is still 

an open question for which value (s) of this parameter the invariants ()5.10|) can be 
'trigonometrized' leading to an algebraic form of the trigonometric E 6 model. 

Finally, the rational Hamiltonian h^yi . . . y§) can be written as 

ftgfel ..., 6)=A , t _^_ + B ^, (5.14) 



where summation over a, b = 2, 5, 6, 8, 9, 12 with a < b is carried out with the 
coefficient functions: 

^2,2 = 8r 2 , A 2 , 5 = 20r 5 , A 2fi = 24r 6 , 

A 2 ,s = 32r 8 , .4.2,9 = 36r 9 , ^2,12 = 48ri 2 , 
A 5 , 5 = 4r 8 , ^5,6 = 54r 9 + 243r 2 2 r 5 , A 5 , 8 = 48r 5 r 6 + 90r 2 r 9 , 
A 5>9 = 12ti 2 + 162r 2 r 5 2 , A 5> i 2 = 36r 6 r 9 + 81r2r 5 r 8 , 
A,6 = 1080r 5 2 + 270r 2 r 8 + 162r 2 2 r 6 , A 6)8 = 144r_ + 2754r 2 r 2 + 324r|r 8 , 
A,9 = 234r 5 r 8 + 405r 2 2 r 9 , A 6)12 = 567r 2 r 5 r 9 + 72r 8 2 + 540r 5 2 r 6 + 486r 2 2 r 12 , 
„4 8j8 = 4374r 2 r 2 + 48r 6 r 8 + 504r 5 r 9 + 216t 2 ti 2 , 
A 8 ,9 = 540r 5 3 + 72r 6 r 9 + 378r 2 r 5 r 8 , 
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A,i2 = 729r 2 2 r 5 r 9 + 162r 9 2 + 72r 6 r 12 + 270r 5 2 r 8 + 108r 2 r 2 + 972r 2 r 2 r 6 , 
„4 9)9 = 216r 2 r 5 r 9 + 12r 8 2 + 144r 2 r 6 , 
^9,12 = -162r 2 r 5 r 12 + 270r 2 r 9 + 189r 2 r 8 r 9 + 144r 5 r 6 r 8 , 

^12,12 = ^r|r 9 2 - 324r 2 r 5 r 6 r 9 + 36r 6 r 8 2 - 648r 2 r 12 + 432r 2 r 2 

+594r 5 r 8 r 9 + 162r 2 r 8 r 12 , (5.15) 

and 

B 2 = -Aujt 2 + 24(1 + Ylv) , B 5 = -10wr 5 , B G = -12cur 6 + 405(1 + 6z/)r 2 2 , 
23 8 = -16ujt 8 + 96(1 + 3z/)r 6 , £ 9 = -18wr 9 + 216(2 + 3z/)r 2 r 5 , 

B X2 = -2Aujt 12 + 108(7 - 6z/)r 2 + 324(1 + 2z/)r 2 r 8 . (5.16) 

There is a one-parametric algebra of differential operators (in six variables) for 
which p^ 1 ' 1 ' 2 ' 2 ' 2 ' 3 ) ( see fl^TT)) 

is a finite-dimensional irreducible representation space. 
Furthermore, the finite-dimensional representation spaces appear for different inte- 
ger values of the algebra parameter. They form an infinite non-classical flag which 
coincides with V^ E ^ (j5.11|) . We call this algebra e^. Like the algebras g^ 2 \ p 4 > 
introduced in 0) and jH], respectively, in relation to the G 2 and F 4 models, the 
algebra e*- 6 -* is infinite-dimensional yet finitely-generated. It will be described and 
studied elsewhere. The rational Eq Hamiltonian in the algebraic form ()5.14j) with 
coefficients (I5.15J) . 1)5.16)1 can be rewritten in terms of the generators of this algebra. 

The operator (j5.14j) is triangular in the basis of monomials t 2 Vf VfVgVfl. 
One can find the spectrum of ()5.14j) . h!^(p = —2eip, explicitly 

eni,n 2 ,n 3 ,n4,n 5 ,n 6 = uj(2th + 5n 2 + 6n 3 + 8n 4 + 9n 5 + 12n 6 ) , (5.17) 

where n« are non-negative integers. Degeneracy of the spectrum is related to the 
number of partitions of an integer number n, n = 0, 1, 2, ... to 2n\ + 5n 2 + 6n 3 + 
8n 4 + 9n 5 + 12n 6 . The spectrum does not depend on the coupling constant g, it 
is equidistant and corresponds to the spectrum of a set of the harmonic oscillators. 
Finally, the energies of the original rational E G Hamiltonian ()5.1|) are E = E + e. 
As an illustration the first eigenfunctions are presented in the Appendix B. It is 
worth noting that the Hamiltonian (|5.14|) possesses a remarkable property: there 
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exists a family of eigenfunctions which depend on the single variable r 2 . These 
eigenfunctions are the associated Laguerre polynomials. This property admits to 
construct a quasi-exactly-solvable generalization of the rational Eq model, as will be 
done elsewhere. Due to enormous technical difficulties we were unable to describe 
explicitly the boundary of the configuration space in the Weyl- invariant variables 
r's similarly to what was done in G 2 and F 4 cases. 



6 The rational Ef model 

The Hamiltonian of the rational E 7 model is built using the root system of the 
exceptional E 7 algebra. A convenient way to represent the Hamiltonian is to write 
it in the 8— dimensional space {x±, x 2 , ■ ■ ■ x$} and impose the constraint x$ = —x 7 , 



(6.1) 



where a; is a frequency and the root generated part of the potential depends on a 
single constant g = v{y — 1): 



Ve 7 =9 

j<i=l 



+ 



+ 9 



1 



(x 7 - x 8 y 



t [l (-xs + xr-^i-lpxjj 



2 ' 



(6.2) 



with Uj — 0, 1, and $^j=i Vj = odd. The configuration space is given by the principal 
E 7 Weyl chamber. 

Let us introduce the new variables 



Vi 



Xi , -i = 1 ... 6 



(using the constraint y 7 = 2x 7 ) 



y 7 = x 7 - X 8 , 

Y = -{x 7 + x 8 ) , (using the constraint Y — 0) 



In these variables the Laplacian becomes 



ft 2 1 f) 2 

A (8) _ A (6) , , 

~ y + d y y 2dY 2 ' 



(6.3) 
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and the potential part of (|6.1j) depends only on {y% . . . y 7 }: 
V 



10 

T 



i=l 
6 



1 1 

+ 



1 [(yi + yj) 2 (yi-yj) 2 \ vi 



i 



^ [J (yr - E- =1 (-i; 



2 ' 



(6.4) 



In this formalism imposing constraints means the restriction to the eigenfunctions 
having no dependence on Y. Hence, F-dependent part of the Laplacian, the last 
term in ()6.3j) . can be dropped off. 

The ground state eigenfunction has a form 



where 



* = (A^A^yrYA^e 



-~^' : -" if+ '' I < E = 7 -uj(1 + 18u) 



(6.5) 



,(6) 



I [ (yi ± %) 



j<i=i 



^ 7 = n u+ec- 1 )"^ 

with z/j = 0, 1 and $^ =1 = odd and g = — 1). 

In order to find variables leading to algebraic form of gauge-rotated Hamiltonian, 



h { £(yi...y 7 ) = -2y - 1 (HE?-E )(y 1 ...y 7 )y 



(6.6) 



let us take the Weyl-invariant polynomials, obtained by averaging over the 56- 
dimensional orbit Q generated by the vector (e 7 — e 6 ), 



56 



(6.7) 



k=\ 



(cf. (j2.7|) ). where a = 2, 6, 8, 10, 12, 14, 18 are the degrees of the E 7 invariants and 
ujk-, k = 1, 2, ... 56 are the orbit elements. The orbit variables ta are functions of 
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Ui . . . y-j. The invariants of a fixed degree are defined ambiguously, up to non-linear 
transformations, similar to (|3.8j) . (|4.7j) . (|5.9j) 











- 4 Q) +A( 6 )(4 n) ) 3 , 




4 n) - 


-> 4 Q) +4 8) 4 n) 4 n) +4 8) (4 n) ) 4 , 




-(H) 

no H 


- 4? + 4 10) (4 n) ) 2 4 n) + 4 10) 4 c) 4 n) + 4 10) (4 n) ) 5 




^12 ^ 


- 4? + A^it^f + A^tfh^ + Af } (4 n) ) 2 4 n) 

+4 12) (4 n) ) 6 , 


+4 12) (4 n) 



(6.8) 



6 



(Q) , (fi) 4 (14).(0)/.(n)N 2 , 4 (14).(0).(0) .(14)/.(n)^2.(n) ^i*,/ jimj, 

^14 l— r 14 "+" 71 1 r 2 l r 6 J "+" A 2 l 2 l l2 "+* ^3 l r 2 J r 10 + ^4 l r 2 J 1 
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+4 14) 4 n) 4 n) + 4 14) (4 Q) ) 4 4 n) + 4 14) (4 n ^ 7 



.(«) .(o) , 4 (is)/.(n)N 3 , 4 (i8).(n).(n) , 4 (i8).(n).(n) , ^iB;.yi;/. 

r 18 l— *■ r 18 + A l l r 6 J + ^2 r 6 l 12 + ^3 r 8 r 10 + ^4 r 2 l r 8 7 

(is) .(n) ,(fi) ,(fi) , ,(i8)/.(fih2,(fi) , 4 (i8)/.(n)v2 + (n).(n) , 4 (is)/.(n)x3 + (n) 



+^5 r 2 r 6 r 10 "+" ^6 l r 2 J r 14 + ^7 l r 2 J r 6 r 8 + ^8 l r 2 ) l 12 

'-6 



+4f(4 n) ) 9 • 

Our goal is to tune the parameters A's so as to get the algebraic form of the Hamilto- 
nian (if it exists) and a minimal flag. After very cumbersome analysis we discovered 
a two-parametric set of variables (for simplicity we omit subscript (Q) in variables 

4 n) ), 

1 16 1 3 16 16 1 4 

To — — to , Tr — te. -\- tn , Tq — to -)- totp. — tn , 

3 3 108 2 ' 5 45 2592 2 ' 

_ 64 4 1 5 1 2 

Tin — tin — tote, — to ~f~ totfi . 

u 315 105 466560 2 405 2 

_1024 5 6 64 3 184 2 (12) 256 2 

Tl2 "^5~ tl2 + 419904 t2 ~ 3645 Me + 405* 2 * 8 " ^ Ml ° " 405* 6 ' 
_4096 43 4 832 2 256 202 s 

T14 — tiA — totp. -\- totp. — tfita — tota 

2233 634230 2 35235 6 1305 246645 2 

341 7 18176 8768 2 

5114430720 2 43065 246645 2 ' 

_ 262144 635120768 2 534492928 3 2 49527971 6 

Tie — — — tie ~\- —~ totfitc — tot fi -(- tfito 

3687 505211175 2 5846015025 2 6 210456540900 2 

4002704 - 192754688 3 9332528 , 2735503 q 

5846015025 2 8 285805179 2 12 233840601 2 10 16971215458176 2 
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262144 3 54099968 (18) 38912 2 

+ 1493235* 6 ~ 17421075 Wl ° " ^ ~ 30725^^ 

458752 98099200 , , . 

*6*12 ^9*14 , (6.9) 

55305 24699213 2 ' V ; 

(12) (18) 

where A 2 , A% are parameters, leading to an algebraic form of the Hamiltonian 
and in which the flag is minimal. We denote the minimal flag V^ E7 \ This flag is 
generated by the polynomials 
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^(1,2,2,2,3,3,4) _ / P2 pe J8J10J12 1 n < r, _|_ 9n _|_ 9n 

' n — \ T 2 r 6 T 8 T 10 T 12 T 14 T 18 I U ^ T>2 + ^6 + 

+ 2 Pl0 + 3p l2 + 3 Pl4 + Ap ls < n) , (6.10) 

with the characteristic vector 

/ = (1,2,2,2,3,3,4) , (6.11) 

fcf. fl3.13j) . ()4.9|) . ()5.12|) ). It again coincides with the highest root of the E 7 algebra 
confirming the Kac's conjecture. Therefore Vn = p^ 1 ' 2 ' 2 ' 2,3 ' 3,4 ). 

It is worth to emphasize that the most general polynomial transformation which 
preserve each linear space T^n ^ is 

sa — > s 6 + a 6 sl + b 6:1 s s + 6 6i2 si , 

S S — *• S S + a§S2 + &8,l s 6 + bs^SlO j 
s 10 ~~ * s 10 + ^lO- 5 ^ + ^10,l s 6 + &10,2-S8 , 
«12 — > «12 + a 12 s 2 + bl2,lS 2 S 6 + &l 2 , 2 S 2 S 8 + b 12 ^S 2 S 10 + C 12 S U , 
S M — > S14 + a 14 S2 + &14,lS 2 S 6 + &14,2S 2 S 8 + &14,3-S2-SlO + C 14 S 12 , 
Sl8 " > Sl8 + disst + bl8,isls 6 + b 18 , 2 sls S + &l 8 ,3-S 2 -Si + Ci 8i iS 2 Sl2 + Ci 8>2 S 2 Si 4 

+^18,l s 6 + ^18,2-Sg + ^18,3 S 10 + ^18,4-S6 S 8 + ^18,4-S6 S 10 + ^lS^S-SlO j (6-12) 

where {a, b, c, d} are arbitrary numbers. It is worth to mention that there exist two 
exceptional sub-transformations which are common for (|6.8|) and (J6.12|) - when we 
vary the parameter 6 123 in s 12 (A 2 12 ^ in t^) and dis^ in s 18 (A^ in t^$) keeping 
all other parameters fixed. Thus, we conclude that there exist a two-parametric set 
of invariants of the fixed degrees leading to algebraic form of the Hamiltonian 

and moreover preserving the minimal flag. Hence the parameters and A^ 

can be chosen at our will. We set them equal to zero, A { 2 12) = 4 18) = 0, which 
fixes the coefficient functions in the algebraic form of the Hamiltonian (see below 
(j6.13|0 in the form of polynomials of lowest degrees. Since an algebraic form of 
the trigonometric Ej model is not known so far, it is an open question for which 
value(s) of these parameters the invariants ()6.8|) would be 'trigonometrized' leading 
to an algebraic form of the trigonometric E 7 model (if such a form exists). 



28 



Finally, the Hamiltonian can be written as 

h% = A ab7 ^— + B a ^- , (6.13) 

7 OT a dr b dr a 

where a, b = 2, 6, 8, 10, 12, 14, 18 and a < b with the coefficient functions: 
A 2 , 2 = 16r 2 , A 2 , 6 = 48r 6 , A 2 ,s = 64r 8 , A 2 ,io = 80r 10 , 
A12 = 96r 12 , .42,14 = H2ri4 , A,i 8 = 144r 18 , 

A,6 = -8r 2 2 r 6 - 160r 2 r 8 + 17280r 10 , A 6 , 8 = 1920r 2 r 10 + 96r 12 - 16r 2 2 r 8 , 
A,io = -H2r 14 - 24r 2 r 10 , A 6 , 12 = -128r 8 2 + 3456r 10 r 6 - 1440r 2 r 14 - 24r 2 2 r 12 , 
A,i4 = -128r 8 r 10 - 36r 18 - 32r 2 2 r 14 , A,is = -23040r 2 r 2 - 40r 2 r 18 

+ 3200r 8 ri4 - 4224t 10 ti 2 , 
A 8 ,8 = 336ri4 + 288r 2 2 rio + 8r 6 r 8 + 12r 2 Ti 2 , A 8 ,io = -16r 2 Ti 4 + 16r 6 r 10 , 
A 8 ,i2 = -1344r 8 r 10 + 480r 2 r 6 r 10 - 108r 18 - 16r 2 r 2 + 12r 6 r 12 - 192r 2 2 r 14 , 
A,i4 = -960r 2 - 16r 2 r 8 Ti + 20t 6 ti 4 - 5t 2 t 18 , 

A,i8 = -3456r 2 2 r 2 - 576t 2 t 10 t 12 + 448r 2 r 8 ri 4 + 24r 6 Ti 8 - 9984ti Ti 4 , 

4t 8 Tio Ti 8 2 

Ao,io = — 7, 1- ~r > -4-10,12 = 960r 10 - 32r 2 r 8 Ti - 20r 6 ri 4 + t 2 Ti 8 , 

3 6 

Ao,i4 = ~ T g Tl4 + 4ri Ti 2 , Ao,i8 = 288r 2 ri ri4 - 871 ^~ 18 + 1Qt 12 t u - 384r 6 r 2 , 

^4.12,12 = H520r 2 r 2 + 192r 2 ri r 8 - 24t 2 t 8 Ti 2 + 36r 2 2 ri 8 - 336t 2 t 6 Ti 4 

- 16r 6 r 8 2 - 1056r 8 r 14 + 1152r 10 r 12 + 576r 6 2 r 10 , 

»4i2,i4 = 288r 2 2 r 2 + 24r 2 ri ri 2 - 64r 2 r 8 r M - 16t 6 t 8 t 10 - 6t 6 t 18 + 1824ri ri 4 , 
Ai 2 ,i 8 = 2688r 2 4 + 73728r 2 2 r 10 ri 4 - Q72t 6 t w t 12 + 544r 6 r 8 r 14 + 384r 2 r 12 r 14 

- 288r 2 r 8 r 18 + 9216r 8 r 2 + 2592r 10 T 18 - 27648r 2 r 6 r 2 , 

4.14,14 = 32r 2 Ti Ti4 - T g Tl8 + 4ri 2 Ti4 - 64r 6 r 2 , 

A 14 ,i 8 = 23040-4 + 64r 2 ri4 + 128t 2 t 2 4 - 384r 2 r 2 r 8 - 608ti t 6 ti 4 

+ 56t 2 ti ti 8 - 64r 8 ri ri 2 + 4ti 2 Ti 8 , 
A 18 , 18 = 55296r 2 2 r 1 3 - 768r 6 r 10 r 18 + 64r 8 2 n 8 + 285696t 2 ti 4 + 4608r 6 r 8 r 2 

+ 11520r 2 r 2 ri 2 + 384ti t 2 2 - 640r 6 r 2 4 + 128t 2 Ti 4 Ti 8 - 640r 8 ri 2 ri 4 

- 15360r 2 r 8 ri ri4 , (6.14) 

and 

B 2 = - Aujt 2 + 56(1 + 18z/) , B 6 = -12lut 6 - 24(1 + 10z/)r 2 , 
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16ut 8 + 20(1 + 6z/)r 6 , B 10 = -20ujt 10 - 2(1 + 2z/)r 8 , 
24cjt 12 + 480(11 + 18i/)ti - 16(5 + 18z/)r 2 r 8 , 
28cjt 14 + 6(1 + 2z/)r 12 + 24(3 + 2u)t 2 t w , 

36cur 18 - 32(49 + 78z/)r 10 r 6 + 32(29 + 18z/)r 2 r 14 + 64(1 + 2z/)r 8 2 . 

(6.15) 

There exists one-parametric algebra of differential operators (in seven variables) 
for which pp< 2 > 2 < 3 > 3 ' 4 ) ( se e (JfTTUJ)) is a finite-dimensional irreducible representation. 
Furthermore, the finite-dimensional representation spaces appear for different inte- 
ger values of the algebra parameter. They form an infinite non-classical flag which 
coincides with V^ 7 ^ ()6.10|) . We call this algebra e^ 7 \ Like the algebras g^ 2 \ p 4 ' 
introduced in jB]) and jS], respectively, in relation to the G 2 and F 4 models, the 
algebra e^ 7 ' is infinite-dimensional yet finitely-generated. It will be described else- 
where. The rational E 7 Hamiltonian in the algebraic form ()6.13|) with coefficients 
f|6.14|) , ()6.15|) can be rewritten in terms of the generators of this algebra. 

The operator (j6.13J) is triangular in the basis of monomials rfVf VfVfoVf 2 1 Vf 4 14 Tf 8 1 
One can find the spectrum of (J6.13j) . (p = —2e(p, explicitly, 

eni,n 2 ,n3,n4,n 5 ,n 6 ,n 7 = 2 u (n x + 3n 2 + An 3 + 5n 4 + 6n 5 + 7n 6 + 9p 7 ) , (6.16) 

where are non-negative integers. Degeneracy of the spectrum is related to the 
number of partitions of integer number n, n — 0, 1, 2, . . . to ni + 3n 2 + 4n 3 + 5n 4 + 
6115 + 7tiq + 9n 7 . The spectrum does not depend on the coupling constant g, it 
is equidistant and corresponds to the spectrum of a set of harmonic oscillators. 
Finally, the energies of the original rational E 7 Hamiltonian (J6.1|) are E = E + 
e. As an illustration the first eigenfunctions are presented in the Appendix C. 
The Hamiltonian ()6.13|) possesses a remarkable property: there exists a family of 
eigenfunctions which depend on the single variable r 2 . These eigenfunctions are 
the associated Laguerre polynomials. This property allows to construct a quasi- 
exactly-solvable generalization of the rational E 7 model. It will be done elsewhere. 
Due to technical difficulties we were unable to find explicitly the boundary of the 
configuration space (in other words, the boundaries of the E 7 Weyl chamber) in the 
Weyl-invariant variables r's similar to what was previously done for G 2 and F 4 cases 
(see Section 3 and 4, correspondingly). 



B n = - 
B u = - 
B 18 = - 



30 



7 The rational E% model 



The Hamiltonian of the rational E 8 model is built using the root system of the 
exceptional algebra Eg. The Hamiltonian is defined in the 8— dimensional space 
{xx,x 2 , ...x 8 } 



H 



Eg 



(7.1) 



i=l 



where u is the oscillator parameter and Ve 8 is the root-generated potential with the 
coupling constant g = v{y — 1): 



V Es = 9 Yl 



j<i=i 



+ 



_(xi -\- Xj) (Xj Xj) 



[l {** + Ej=i(-D"^ 



2 ' 



(7.2) 



where Uj = 0, 1 and Ylj=i u j = even - The configuration space is given by the 
principal E 8 Weyl chamber. 

The ground state eigenfunction of the Hamiltonian is of the form 

= (Af A^YA Ea e-^=^ , E = + 30u) , 

where 



(7.3) 



.(8) 



I I {Xi ± 



X 



j<i=l 



3 1 ' 
7 



Wj} V 3=1 

with Vj = 0, 1 , X/j=i = even an d g = v{y — 1). 

In order to find variables leading to algebraic form of gauge-rotated Hamiltonian, 

h^{ Xl ...x 8 ) = -2^ - 1 {He 8 -E ){x 1 ...xs)^o , (7.4) 

let us define a basis in the form of Weyl-invariant polynomials, averaged over one of 
the smallest orbit, of length 240, generated by some positive root, 



240 



u k ■ x 



ujk G 0(a positive root) , 



(7.5) 



k=l 
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(cf. f)2.7p ). where a = 2, 8, 12, 14, 18, 20, 24, 30 are degrees of the lowest E 8 invariants 
and ujk, k = 1, 2, . . . 240 are the orbit elements. The orbit variables 4 are functions 
of X\ . . .x%. In general, the invariants of fixed degree are defined ambiguously, up to 
a certain non-linear transformation, cf. ()3.8|) . ()4.7j) . ()5.9|) . ()6.8|) 

fc 2 — ' L 2 ) 

t w _> 4 n )+A( 8 )(4 n) ) 4 , 

.(o) .(O) 4 (i2)/ ,(nu 2 (n) ,(i2W,(n)s 6 

r 12 ~~ r 12 ' A l \ l 2 ) l 8 ' 2 l r 2 J ) 

(Q) .(14) (Q) (Q) 4 (14W,(n)s 3 ,(n) , 4 (14)/,(n)^ 7 
^14 ~~ r l4 + A l l 2 l 12 ' A 2 \ l 2 ) l 8 + ^3 l r 2 



.(a) (fi) 4 (is) , (n) / , (n) \ 2 , 4 (i8)/,(n)^2.(n) 4 (i8)/,(nK 3 ,(n) 4 (i8)/.(n)x 5 .(n) 

l 18 ~ * l l8 "+" ^1 r 2 l r 8 J "+* ^2 l r 2 J r 14 + ^3 \ l 2 ) l l2 "+* ^4 l r 2 J r 8 

+4 18) (4 

.(o) .(n) 4 (20),(n) .(o) 4 (20),(n),(n) 4 (20),,(n)s 2 ,,(n)x2 , 4 (2o),,(n)^ 3+ (n) 

r 20 ~~ * l 20 ' A l l 8 l 12 + ^2 r 2 r 18 + ^3 \ l 2 ) \ l 8 ) + ^4 l r 
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+Af(4" ; )M? + 4 20) (4"TC + 4 ZUj (* 



2 



,(n) (n) 4 (24)/.(n)x2 , 4 (24)/.(0)n 3 , 4 (24).(n).(n).(n) , 4 (24),,(fi)s 2 ,(n) 

t-24 ~~ * t 24 "r ^1 ^12 J ^2 l 6 8 J ^3 L 2 L 8 L IA ^4 14 J t 20 

+At i \tf ) ftfh { s + 4 24) (4 Q W + 4 24) (4 n) ) 4 (4 n) ) 2 + 4 24) (4 n) ) 5 4? 
+4 24) (4 n) ) 6 *£ ) + Ag 4) (4 Q) ) 8 4 n) + A^\tfY , 

AO) _ ,(n) 4 (30).(o).(n) 4 (30) r (fi)x 2 .(n) 4 (30) .(«)/. (fl)v 2 , 4 (30).(n),(n).(n) 

^30 ^30 ^1 t 12 18 "r ^2 ^8 I "14 i" ^3 t 2 ^14 J "r ^4 t 2 6 8 t 20 

4 (30),(n),(n) 2 ,(n) 4 (30)/,(n)x 2 .(fi).(n) 4 (30),,(n)N 2 ,(fi) # (n) , 4 (3own) v3 (n) 

"l - ^ t 2 ^8 fc 12 ^6 ^2 J t 12 L 14 "I" ^7 l t 2 J L 8 18 ^8 l t 2 J fc 24 

+4 30) (4 n) ) 3 (tS? ) ) 2 + 4f (4 n) ) 3 (4 n) ) 3 + 4f \4 n) ) 4 4 n M? + 4 3 2 0) (*f } ) 5 4? 
+4l 0) (4 n) ) 5 4 n) 42 ) + v4[ 30) (4 n) ) 6 tSg ) + v4[ 30) (4 n) ) 7 (4 Q) ) 2 + 4 30) (4 n) ) s t (n) 
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+Ag 0) (4 n) ) 9 tS? + 4 3 8 0) (4 n) ) n 4 n) + 4 3 9 0) (4 n) ) 15 . (7.6) 

The transformation (|7.6j) depends on 48 parameters. Our goal is to find pa- 
rameters A's such that (i) the Hamiltonian has the algebraic form, (ii) a set of 
polynomial invariant subspaces forming a flag occurs and (iii) the flag is minimal. 
After extremely tedious and cumbersome analysis we discovered a nine-parametric 
set of variables (see discussion below), for simplicity we omit the subscript (Q) in 
variables 4P\ 

T2 = iV 2 ' 

32 



1 13 



8 ~30 8 16200000 2 ' 

16 (12) 2 373 g 

ri2= 21 tl2 + A Ms+ 2551500000 * 2 ' 

64 2531 7 568 103 3 

7~1 4 — ^14 — ^o — — — totio ~\~ Z ^o ^8 i 

1155 2296350000000 2 51975 1417500 2 

_ 256 (18) 2 3706 3 (18) 2 4051 5 

ri8 "4095 tl8 + A2 Ml4 + 8353125 ^ + ^ *2 f 8 - 1530900000 ^ *s 
330961 9 
+ 8266860000000000 2 ' 



4096 60576512 18641008 



g 

Ton = ton — to tic -|- to t-M 

575025 35320910625 779137734375 2 

1 / 138176 103942624 4 \ 323371 2 2 

+ 38335 V ~ 1575 * 8 ~ 1196015625 * 2 J tl2 + 6538218750 * 2 1§ 
10249681 in 2994007 „ 

76261783500000000000 2 353063812500000 2 8 ' 
(24) 3 /i(24) /t( 24 ) 2 32768 52162303808 3 

T94 = j\o to -|- to to tiA -|- .A,i to ton -)- to4 -)- — — — — — — to tis 

24 2 3 2 8 14 t 4 2 20 t 13101165 ^ 878607651796875 2 

2857817967448 5 197632 2 

~ 3663018665947265625 * 2 14 ~ 362005875 * 12 

1 / 12897642016 2 11785468451047 6 \ 

+ 4367055 V 31134375 * 2<8 + 36777480468750 V * 12 
22206803851 4 2 101271432653 8 

~ 13661608078125000 * 2 * 8 ~ 368863418109375000000 * 2 * 8 
38183226373283 12 

+ 8764194814278750000000000000 * 2 ' 

- _j(30) / 2 / , 4 (30) . ,2 , 4 (30) . . . 4194304 

^30 ~A 2 t 8 tl4 + ^3 J 2 tl4 + A 4 % 2 H t 2 Q + ^4439375 ts0 

87361458176 3 98943157092328832 5 



412072580390625 2 24 14540238544445947265625 2 20 

1 / 82608128 76042100276224 2 44731593575760671656 fi 

7632625 V ' 945 12 14671762875 2 8 6158701713076171875 2 J ° 18 
1 / 387292030976 2 

+ 7632625 V 128638125 2 12 

8358520483853754832 4 621999713517306328312 „ , , 

+ - n1 , », omonnTO10 r ^2 *8 + otrtro^oonnooi -^^o-^ ^2 J *14 



101437439980078125 z 8558783998319091796875 
162900154624 3 2 1 /29318297174541 

81160947912515625 * 2 tl2 + 7632625 V 271859135625 
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79803090068621091328 5 101275905787214796443 9 \ 

4564684799103515625 2 8 ~ 15483227333642578125000 V * 12 

10147539336141079 , , 400141317989263534193 7 » 

2 8 i7c;/onoQconnt:^oinc/fi07cnnnnnn 2 8 



1082844671787333984375 z 17542083682954810546875000000 
105630144697706078621 n 

+ 33831 161388555706054687500000000 2 8 

7373766632847391460197357 15 

~ 146296767815759210806406250000000000000000 * 2 ' ^ ' ' 

which lead to the flag which we think is minimal. We denote this flag V^ Ea \ The 
flag p( £s ) is generated by the spaces of polynomials 

^(1,3,5,5,7,7,9,11) = {rftf^ffTtfTtf-ffT*"] < p 2 + 3p 8 + 5 Pl2 + 5 Pl4 + 7 Pl8 

+ 7p 20 + 9p 24 + llpao < 7i> , (7.8) 

with the characteristic vector 

/ = (1,3,5,5,7,7,9,11), (7.9) 

(cf- dmD , (BSD , flnH) , ijnrrt ). Hence = p(i>3,5,5,7,7,9,u)_ The characteristic 

vector does not coincide with the highest root /highestroot = (2,2,3,3,4,4,5,6) sug- 
gested by Kac. We were not able to find variables in which the flag /highestroot would 
be preserved by the rational E$ Hamiltonian. 

The most general polynomial transformation which preserves each linear space 
Pf'™ 11 ',^ 0,1,2,... is 

S2 — > S 2 , 
s 8 ~~ * a 8,l s 8 + a 8,2S2 > 

Si 2 — > ai2,lSl2 + 012,2^14 + dU^Ss + 012,4^3 , 

2 5 
S14 — > Ou^Su + 014,2^12 + Ol4,3S2 s 8 a 14,4S 2 , 

1 1 2 1 2 1 2 1 4 1 7 

Sis — > ai8,lSl8 + Oig 5 2S20 T a 18,3 s 2'S 8 + «18,4 S 2 S 14 T Ol8,5 s 2 S 12 + a 18,6 s 2 S 8 + a 18,7 s 2 ) 

2 2 2 4 7 

S20 O20,l s 20 + ^20,2^18 + CL20,3S2S 8 + a20,4<S 2 Sl4 + O20,5S 2 Sl2 + a20,6S2 S 8 + ^20,7*2 ' 

3 2 2 

S24 — ► 024,lS24 + 024,2<S 8 + 024,3S2S8Sl4 + Cf24,4>S2S8Sl2 + a24,5"S 2 S20 + a 24,6"S 2 Sl8 

+a 2 4,7 s 2 S 8 + a 24,8'52'5l4 + C-24,9'S2'5l2 + <224,10 S 2 S 8 a 24,llS2 , 
s 30 — y °30,l s 30 + a 30,2 s 8 s 14 + a 30,3 s 8' S 12 + a 30,4S2Si4 + O30,5 s 2 s 12Sl4 + Cl30,6 s 2 s l2 
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2 2 3 3 

+ a 30,7 s 2-58'S20 + Q30,8 S 2 S 8 S 18 + a 30,9 s 2 S 24 + 030,10-32-58 + &30,11 S 2 S 8 S 14 

+ a 30,12-52- s 8 s 12 + C t 30,13- S 2' s 20 + O-30,U S i S 18 + a 30,15 s 2 S 8 + a 30,16 s 2 s 14 

+030,17^2512 + a 3 o,i8S2S 8 + a^igsl 1 ,(7.10) 

where a's are parameters. There exist nine exceptional sub-transformations for 
which ()7.6|) and ()7.10|) coincide - when we vary the parameter 012,3 i n $12 (^i 12 "* in 
4^) and ai8 3 and ai 8 4 in s 18 (A^ 18 ^ and A^ 18 -* in 4?)> correspondingly, and a 2 4 2, a 24 3 
and a 24 , 5 m s 24 {Af\A^ and 4 24) in 4?), correspondingly and a 30 ,2,a 3 o,4 and 
a 30 ,7 in s 30 (^4 2 3 °\ ^4 3 °^ and in 4o' ) ); correspondingly Of course, this coincidence 
appears when all other parameters of the transformation are kept fixed. Thus, 
one can draw a conclusion about existence of nine-parametric set of invariants of 

(r) 

the fixed degrees leading to algebraic form of the Hamiltonian h E& and moreover 
preserving the flag p( 1 ' 3 ' 5 ' 5 ' 7 ' 7 ' 9 - 11 ). Hence, the parameters can be chosen by following 
our convenience. In a simply-minded way we set all of them equal to zero, 

.(12) _ ,(18) _ .(18) _ ,(24) _ .(24) _ ,(24) _ ,(30) _ ,(30) _ ,(30) _ n 
A l — A l — ^-2 — ^-2 ~ ^3 — ^3 — ^2 ~ ^3 — ^4 ~~ u ) 

which fixes the coefficient functions in the algebraic form of the Hamiltonian (see 
below (J7.11j0 in form of polynomials of lowest degrees. It is an open question for 
what value (s) of these parameters the invariants (|7.6j) are 'trigonometrized' leading 
to an algebraic form of the trigonometric E$ model, if such a form exists. 

Finally, the Hamiltonian h^\x\ . . . x§) can be written as 

q2 q 

hP(x 1 ...x 8 )=A ab —— + B a — , (7.11) 

8 dT a dn dr b 

where a, b = 2, 8, 12, 14, 18, 20, 24, 30 and a < b with the coefficient functions: 
A 2 ,2 = 8r 2 , A 2 ,s = 32r 8 , A 2 ,u = 48ri 2 , A 2 ,u = 56r M , 

A.is = 72ri 8 , A 2 ,2o = 80r 20 , -4 2 ,24 = 96r 24 , -4 2 ,30 = 120r 30 , 

21 7 , 27 5 2421 2 14664 2 12 3 
-4 8 ,s = -^t u +—t 2 t 12 , -4 8 ,i2 = — — t- 2 t 8 + ^— r 2 Ti4 + 81ri 8 — r 2 r 8 -— r 2 ri 2 , 

32 2 , 2349 2 16640 3 10 5 

-4 8 ,i4 = 75r 20 + Y^ r s Tl2 + g r 2 r i8 , -4 8 ,i8 = — — r 2 r 20 + r 8 - 18t 2 t 14 , 

3744 , 2 88 2 16 2 12232 27 , 

+270r 24 + ^- r 2 r s + ^ r 2 r 8 r i2 + ^ r i2 ^T- 2 r 8 ri4 - —t$t 1b , 
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3 32 . 43264 3 3 2288 4 16 2 

A» = 5^r 24 + ^r 12 r 14 , A 8 , 2 , = r 30 - ~^r 2 r 8 + + ^t 2 t 12 t u 

832r 2 2 r 8 r 18 64r 12 r 18 126r|r 20 39r|r 24 9936r 2 r 8 r 20 311168t 8 2 t 1 4 352r 2 r 1 2 4 
+ 1575 + 1125 5 5 35 + 19845 + 125 ' 

. 48 1555840 2 27 3 43264 2 43264 2 

A,30 = yr 12 r 24 - ^^r 8 r 20 - -r 2 r 30 + ^r 8 r u + ^ 
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22275 1185120 3 2 810 7 
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2688 innnn 4415680 2 , 2218185 4 rrtJ 2 

^12.18 = — T 2 Ti 2 Ti 4 + 12600r 8 T 20 — T 2 T § — T 2 T 20 + 504r 14 

5 49 14 

270 7 72 2 2 53640 6 2 264 4 1891584 3 
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10112 2 20040 81 5 381375 2 
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7 7 7 z 14 

48 4 144 2 9 3 4800 315008 
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384 2 
+ ^^ 2 r 14 , 

rw , no 10112 3 2 43264 5 3 1024 2 117 5 

^ 1224 = 54r 2 7 r 20 + 408r 14 r 20 + —^r-T»Tf 4 + -^r^s + ~^~n fr 18 + — r|r 24 

1504 4 21408 3 _ 11152 32 4 64 6 

— ^7^T 2 T8^i8H 7^— T 2 T 8 r 20 -576r 2 ri 2 r 2 o+ t 2 ti 4 ti 8 + — — t 2 ti 2 ti 4 - — r 2 r 8 ri 4 

too 4y 1/0 0Z0 ( 

60764416 2 2 3328 96 2 6600 1968 2 

"^6305- r27 " 8ri4 + T05" T8ri2Tl4 " ITS" 2 " 12 " 18 + ~ —™ ' 

3056 4 36000 6 _ 2 1730560 4 4 131456 5 2 
^12,30 = — T- 2 ri 2 r 20 ^— r 2 r 8 r 20 + 27000r 20 926I ^ + 1029 T2T87 " 14 

346112r 2 2 r 8 3 r 12 43264r 2 3 r 2 r 18 9120r 2 4 r 8 r 24 648 2 71552 3 

— r 2 ri 2 r 24 + T 2 r 8 ri 2 ri 4 



3087 9261 49 7 2 945 

31800 81 5 55377920 5 16640 

-t 2 t 8 T3o + 1728r 2 ri 4 r 24 + — r 2 r 3 o H — — — r 8 — r 8 ri 2 r 20 



7 u ou 7 ' " 3969 7 
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469779200 2 2 95352 3 28080 7444736 2 2 

+ T-2 r 8 T 20 —r 2 T 14 T 20 — T 2 T 18 T 20 + 2 205 

1646080 2 263644160 3 43264 
+^^r 8 r 24 - g261 r 2 r 8 r 14 + —^ T ur w , 

„ 3 525 99200 3 11360 2 8 208 

^.14,14 = -1575r 2 r 20 —t 2 t 24 ^— t- 2 t 8 H ^— r 2 r 8 ri 4 - -r 12 r u + — r 8 ri 8 , 

75 4421120r 2 3 r| 1408r 2 r|r 12 23680r 2 4 r 8 r 14 80r 2 2 r 12 r 14 2928 2 
^14,18 = yr 3 o+ m + 21 7 + 3 "Wis 

32 n.rn 5 ™ 3 35080 1902848 2 ^ 2 

+— ri 2 ri 8 + 9450r 2 r 20 + 1575r 2 r 24 H — r 2 r 8 r 20 H — - r 8 ri 4 - 704r 2 r 14 , 

lo / lozo 

. oo 16 32 2 2 61696 2 48 
-4i4, 2 o = o2t 8 t 24 — —Ti 2 t 20 + — r 2 r 14 1323 r2T 8 ri4 + 25 Tl4Tl8 ~~ 2X30 ' 



16 512 2 3 2944 2 1856 2 64 4 2 32 

y^ 2 T- 24 H — — T- 8 r 2 o + 18r 2 r 3 o + yyr 8 r 14 + ^^-r 2 ri 4 ri 8 -yr 2 r 14 + — : 



4352 2 2154496 3 2 9328 4 8696 2 
■yyT 2 r 8 Ti 8 + 6615 r 2 r 8 r i4 ^— t 2 t 8 t 20 —t 2 t 8 t 24 

1664 1384448 2 4 4544 



+ ^25" r2T8ri2ri4 " ^969~ r27 " 8 " —™ t ™ > 
1664 4630912 38600 2 2145894400 5 

^14,30 = —r 2 T 8 T 12 T 20 1984 5 T 2 T 8 r 14T-18 —^20 250047 r2r 8~ 8Tl2T30 

26800 2 2 584r 2 r 8 r 30 1459840 2 726968320 2 o 
+ ^^2^18^0 + 1200r 2 r 14 r 24 + r 2 r 8 r 24 + 83349 r 2 r s r i4 

346112 2 4499456 3 2248832 17667200 3 2 

+ "66T5" r8Tl2Tl4 + TT907" T8ri8 + 216T18T24 " ^323~ r8ri4T2 ° " ~^^ T ^ 

83200 .3. 2 , 7nn ,4„ 
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-T^rgr^ + 7200r 2 V 14 r 20 , 



, o 9 19125440 = 3 70400 9 
^Ii8,i8 = -56700r 2 7 r 20 + 1020r 14 r 20 - 960t 3 t 2 4 ^— rJ ?r| + -^-^ fn 8 

ftilrft 5 15504 4 24064 3 2 559980 3 „„ ifT „ n 

-9450r 2 r 24 H — r 2 r 8 ri 8 H — r 2 r 8 r 12 H — r 2 r 8 r 20 + 68r2r 12 r 2 o - 1 760r 2 r 14 ri 8 

4 134880 6 114368 2 2 1088 464 2 

-256t 2 t 12 Ti 4 H — r 2 r 8 r 14 t 2 t 8 t 14 + yyr 8 ri 2 r 14 - — r 2 r 12 ri 8 

2368 2 164100 166400 4 945 2 

+ U)5" Wl2 H 7 T 2 T 8 T 24 + ^2^8 ' 

16 10880 2 3 1408 2 944 2 667904 3 2 

^18,20 = yT"i 2 T- 24 H — T- 8 r 2 o + 18r 2 r3o + yyr 8 r 14 -— t 2 t 14 t 18 + t 2 t 8 t 14 
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352 4 2 240 4 1088 2 1408 AO 2 5032 

^-T 2 T 14 + -j-T 2 T 8 T 2 o ^T 2 T 8 T 24: +-^rT 2 T 8 T 12 Tu + 48T 2 T 12 T 20 ^— T 2 Ti 4 T 20 , 

, 3328 2 7264 136 2 12784 26172032 5 2 

' 24 = 23625 T2T " 12ri4 — T 2 T 14T24 + — T 2 Ti 2 T 24 ^ T 2 Ti 8 T 20 15435 T 2 T 8 T 14 

26850304 3 3048064 2 2 4432 4 45456 6 647296 2 2 

-T 2 TgTi 4 H T"H r 2 T 8 T 20 7^T T 2 r 14Tl8 H 7^T 2 T 8 T 20 + ———T 2 T 8 T 1A 



83349 8 441 ^ 8 u 175 z 7 z 11025 

692224 23 592 4 94592 8808 3 287744 3 2 

■3^r 2 r 8 r 12 - — r 2 r 12 r 20 + -^r 8 r u r 18 - -^ 2 r 14 r 20 + ^^r 2 r 8 r 18 

256 14344 4 1280 13568 3 

T 8 T 12 T 20 H = T 2 T 8 T 2A + -^-—T 2 T 8 T 30 + T 2 T 8 T 12 Ti A 



45 ° ' 7 A ° 21 3675 

4352 40960 
^r 2 r 8 r 12 r 18 + ^ r . .^ ^ 



!"Tr7^ r 2 r 8 r 12T'l8 H r 8 2r 24 7^T T 2 T 18 ~ 108 T 2 T 30 



47071232 4 4 256 2 4512 , 2 rtftn 2 

+ ^778T T ^ + 1125 Tl2T » + + 200T| » ' 

74054656 2 2 743104 5448 4 692224 2 

As ' 30 = -^6305- r2T8Tl4 + ^T r8ri4T24 " -T^ 7 " 87 " 30 + ^7675 r8ri2ri4 

3328 2 30457856 4 737884160 3 24960 
-^r 2 r 12 r 20 + W12 + ^ r 2 r 8 r 20 - — r 8 r 18 r 20 

448640 3 2 260368000 5 2 692224 2 21112832 2 3 

+— ^-t 2 t 8 t 24 + ^ W20 + ^^ 8 t 12 t 18 - -^"Wis 

10771764224 4 3 216 2 326400r 2 6 r 14 r 20 3757312 5 2 ^ non 4 

-r 2 r 8 ri 4 H — — r 2 ri 2 r 30 1 — — r2r 8 r 14 -7920r 2 ri4r 24 



194481 5 ' 7 441 

79740r 2 r 20 r 24 57600 4 128 2 
1536t 2 ti 4 t 3 o — t 2 t 18 t 20 + — t 12 t u t 20 - 1416r 2 Ti 8 r 24 



294080r 2 2 r 8 r 14 r 20 229924864r 2 2 r 8 2 r 12 r 14 8000r 2 r 8 r 12 r 24 3306368r 2 3 r 8 ri 4 ri 8 
+ 21 416745 + 21 + 2205 

422528r 2 3 r 8 ri 2 r 2 o 95637667840r 2 3 r 8 5 248400r 2 3 r 2 2 19793453056r 8 4 r 14 3354880r 2 r 30 
+ 147 + 1750329 7 + 5250987 + 1323 

4096r 2 r 8 r 2 4 32r 8 r 30 692224r 3 r 14 304r 14 r 24 224r 2 2 r 14 r 2 o 64r 18 r 20 
20 ' 20 525 15 59535 25 5 15 

1664 2 48 32 2 16 68864 

.4-20,24 = — T 2 T 8 T 24 + — Ti 8 T 24 + —T 2 T 8 T^ - 7^ T 12 r 30 + 1575 t 8 t 14 t 20 

1792 „ n 2 512 379136 „ 2 6656 2 

+ 5625 Tl4 " 976T2T2 ° + 3375 r2r8Tl4ri8 + "W^ 14 + 118125 
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21458944 2 3 346112 2 10944 2 5344 4 

-^97675- r2T8Tl4 + 893025 T8Tl2Tl4 " T^ 7 " 2 ™ 7 " 24 " 
544 2 13312 3 2 6656 

+ ^- r 2 T 18T20 T 2 T 8 T 20 ~ -^T 2 T%T 12 T20 , 

32 1008 2 327808 2 12193792 2 2 2 3461120 2 3 
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110080 2 ono.n 4 2 227072 3 11318912 692224 3 

" r 8 r 20 + 30240r 2 r 20 — — t 2 t 8 t 14 t 20 — — — T 2 T 8 T 14 T 2 i + t 8 t 24 
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4499456 2 6656 8893693952 4 71360 2 
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+ 27783 r2rsTl4T24 + ^SSig - T2T8ri8r2 ° + ^3^^™° 
147142400 3 1384448 153620480 4 2 
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147 ^ 6615 147 

13312 452727808 2 6656 

H T"2T"8T"18T"30 19845 r2T8T 14 r 20 — r 2T-12T"20T24 



46419987660800 , 97332232192 - nnin 2 384 2 

+ ^o.oooo r 8 r i8 " 3840r 18 r 2 2 - — T 2 4 r 30 , (7.12) 



and 



992436543 8 47258883 

21 

•180 , . _ 9 

~7 
16 



B 2 = -4ut 2 + 32(1 + 30z/) , B 8 = -16ut 8 - — (1 + 18z/)r 2 3 , 
Bi 2 = -24cur 12 - ^(17 + 207z/)r 2 r 8 + |(1 + 18i/)r 2 5 , 



014 = -28wr 14 - (1 + 30z/)r 12 - -(23 + 354i/)r 2 2 r 8 , 
15 7 
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12 4 8 

B 18 = -3Qujt 18 + —(25 + 282z/)r 2 4 r 8 + -(19 + 198z/)r 2 2 r 12 - -(1037 + 7470z/)r 2 r 14 

(DO 

+ 1^(83 + 435^)r 8 2 , 

r i 9 o *y /? o o 

£ 20 = -40^r 20 + —(2 + 39z/)r 2 r 8 2 - —(1 + 18z/)r 2 r 14 - — (1 + 30z/)r 18 , 
od 75 75 

512 256 16 

B 24 = -48^724+—— (23-195i/)r 2 r 8 r 12 +—— (859+3480z/)r 8 r 14 - — (15331+76230z/)r 2 r 20 
4725 1575 105 



(115 + 6678z/)r 2 V 18 - —(191 + 5590i/)r 2 d r| + —(161 + 4338i/)t 2 4 Ti 4 



1575 ' ' " 735 v " ' 75 ' ' "' ' 

B 30 = -60cur 30 - 76544r2ri2ri4 + — (8935+74646z/)r 2 3 r 8 r 14 +^-(3637+49140;/)r 2 r 8 r 18 
m m 14175 g45 v ) 2 8 19845 V ; 

--^-(128693 + 531090z/)r 8 r 20 + ^^(7 - 30z/)r 8 2 r 12 - ^(839 + 2646z/)r 2 2 r 24 



13312 64 

; (2861 + 37362z/)r 2 2 r|-16(151 + 450z/)r 2 V 20 + — — (5837-5670i/)r 2 4 . (7.13) 



27783 v / z s - v , * 4725 

It can be found one-parametric algebra of differential operators (in eight vari- 
ables) for which pa' 3 ' 5 ' 5 ' 7 ' 7 ' 9 ' 11 ) ( see fl£gj)) 

is a finite-dimensional irreducible repre- 
sentation space. Furthermore, the finite-dimensional representation spaces appear 
for different integer values of the algebra parameter. They form an infinite non- 
classical flag which coincides to (17. 8J) . We call this algebra e^ 8 \ Like the 
algebras g@\ introduced in jUj and [S], respectively, in relation to the G 2 and 
F4 models, the algebra e^ 8 -* is infinite-dimensional but finitely-generated. It will be 
described elsewhere. The rational E$ Hamiltonian in the algebraic form (|7.1H) with 
coefficients (J7.12J) , f)7. 13j) can be rewritten in terms of the generators of this algebra. 

The operator ()7.11|) is triangular in the basis of monomials 
rfrfrgVfiVf^rlfVflVIo 30 . One can easily find the spectrum of (J7TQ), 
h^ip = —2e<p, explicitly 

e ni ,n 2 ,n 3 ,n 4 ,n 5 ,n 6 ,n 7 ,n 8 = 2 w(ni + 4n 2 + 6n 3 + 7n 4 + 9n 5 + 10n 6 + 12n 7 + 15n 8 ) , (7.14) 

where pi are non-negative integers. Degeneracy of the spectrum is related to the 
number of partitions of integer number n — 0, 1, 2, ... to rt\ +4n 2 + 6^3 + 7n 4 + 9^5 + 
10n 6 + 12n 7 + 15n 8 . The spectrum does not depend on the coupling constants g, 
it is equidistant and corresponds to the spectrum of a harmonic oscillator. Finally, 
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the energies of the original rational E 8 Hamiltonian (|7.1|) are E = E + e. As an 
illustration the first eigenfunctions are presented in the Appendix D. 

It is worth noting that the Hamiltonian (J7.11j) has a remarkable property: there 
exists a family of eigenfunctions which depend on the single variable r 2 . These 
eigenfunctions are the associated Laguerre polynomials. This property admits to 
construct a quasi-exactly-solvable generalization of the rational Eg model. It will 
be done elsewhere. Due to enormous technical difficulties we were unable to find 
explicitly the boundary of the configuration space (in other words, the boundaries 
of the E s Weyl chamber) in the Weyl-invariant variables r's similar to what was 
previously done for G2 and F 4 cases (see Section 3 and 4, correspondingly). 

8 Conclusion 

We have found in uniform way that all the rational integrable models associated 
with the root systems of exceptional Lie algebras are exactly-solvable, thus contin- 
uing the analysis of the rational (and trigonometric) systems related to the A n and 
BC n root systems as well as their supersymmetric generalizations carried out in 
Ref . jU Uni I2DI • Our method contains two important ingredients: (i) gauging away 
the ground state eigenfunction and (ii) taking specific Weyl-invariant polynomials 
(functions) as variables. After this procedure each Hamiltonian takes an algebraic 
form thus becoming a linear differential operator with polynomial coefficients. Fur- 
thermore, it has infinitely-many invariant subspaces of polynomials, which form a 
minimal flag. The known characteristic vectors of the minimal flags are collected 
in the Table. The meaning of the Weyl-invariant variables (|3.1U|) . (|3.14j) . (j4.11j) . 
(I4.13|) . (|5.1U|) . (|6.9|) . (|7.7|) in which the flag of invariant subspaces is minimal is 
unclear so far. Namely, why in these Weyl-invariant variables the minimal flag is 
preserved. We show that unlike the rational A n and BC n models which all are 
of the hypergeometric type every rational G2, F4, E G j )8 model is special. Each of 
them is characterized by its own hidden infinite-dimensional algebra which deserve 
a separate investigation. It will be done elsewhere. 

Algebraic forms of the rational G2, F^, ^,7,8 models which we found allow to 
construct quasi-exactly-solvable generalizations j2I] of these models similar to those 
of the paper [22] • It is already done in [23] . It concludes an analysis of rational 
models associated with classical (crystallographic) root systems, while the analysis 
of the rational systems related with dihedral (non-crystallographic) root systems 
i7 3 , i?4, ^{m) is still waiting to be done. 
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The presented work complements previous studies where the algebraic and the 
Lie-algebraic forms as well as the corresponding flags were found for the rational 
and trigonometric Olshanetsky-Perelomov Hamiltonians of A — D series (and their 
super symmetric generalizations) jUE], G2 jH] and F4 jH] models. In order to conclude 
a study of the whole set of Olshanetsky-Perelomov integrable systems appearing in 
the Hamiltonian reduction method it is necessary to perform the same analysis for 
remaining Eqj^ integrable trigonometric models. We consider it as a challenging 
task for future. 
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Model 


Rational 


Trigonometric 


A 


(M__i) 






n 


n 


BC n 


n 


n 


G 2 


(1,2) 


(1,2) 




(1,2,2,3) 


(1,2,2,3) 




(1,1,2,2,2,3) 


(1,1,2,2,2,3) 


E 7 


(1,2,2,2,3,3,4) 


? 


E 8 


(1,3,5,5,7,7,9,11) 


? 



Table 1: Characteristic vectors of rational and trigonometric models associated with 
classical (crystallographic) root spaces 
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A First eigenfunctions of the rational F 4 models 

In this Appendix we present explicit expressions for the first eigenfunctions of the rational 
F4 model at n = 0, 1, 2. 



• n = (one eigenstate) 



n = 1 (one eigenstate) 



00 = 1 , 
eo = . 



h = r 2 - -(6^ + 6^+1) , 

UJ 



n = 2 (three eigenstates) 



4 X) = r| - -(4p + 4i/+l)7s + -^(4/i + 4i/ + l)(6/i + 6i/ + l) , 



,(3) 



1 3 

4 2) =t 6 - — (2/i + 4i/ + l)r| + (2/i + 4i/ + l)(4/i + 4z^ + l)r 2 

+ ^3 (2/z + 4i/ + l)(6/i + 6z/ + l)(4/i + 4*/ + 1) , 



6cj 



4 3) = r 8 --(3^ + l)r 6 + -L(3 I / + l)(2 M + 4i/ + l)r 2 2 
1 



3 -(3i/ + l)(2/x + 4i/ + l)(4/x + 4i/ + l)r 2 

+ ^(3^+ l)(2/x + 4z^ + l)(6/i + 6v + l)(4/i + 4i/ + 1) , 
8a; . 



B First eigenfunctions of the rational Eq model 

• n = (one eigenstate) 

00 = 1 , eo = , 
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n = 1 (one eigenstate) 

6 

0i,i = r 2 - -(l + 12i/) , ei,i = 2w, 
01,2 = T5 , ei,2 = 5w . 

n = 2 (six eigenstates) 

02,1 = r| - — (1 + 9^)r 2 + -^(1 + 9i/)(l + 12i/) , e 2 ,i = 4c; , 

UJ 

02,2 = -4r 6 + — r 2 2 (l + 6z.)-^r 2 (l+6^)(l+9z/) + ^(l + 6^)(l + 9^)(l + 12z/) , 

UJ UJ UJ 

£2,2 = 6w , 

g 

02,3 = (T2 (2 + 9i/))r 5 , e 2j3 = 7w , 

02,4 = r 8 - ^ (1 + 3^)r 6 + (1 + 3u) (1 + 6i/)r| - (1 + 3i/) (1 + 6i/) (1 + 9z/)r 2 
9720 

+-^4-(! + ^X 1 + H(! + 9i/)(l + 12i/) , e 2j4 = 8cj , 

54, N 216, 
02,5 = r 9 (2 + 3^)r 2 r 5 + —=-(2 + 3i/)(2 + 9z/)r 5 , e 2 , 5 = 9kj , 

02,6 = rf - -r 8 + — (1 + 3^)r 6 - ^(1 + 3^(1 + 6^)r 2 2 

UJ UJ Z UJ a 

3240 3888 
+-^ 3 -(1+3i/)(1+6i/)(1+9i/)t5 — ^g-(l+3i/)(l+6i/)(l+9i/)(l+12i/) , e 2 , 6 = IOcj . 

First eigenfunctions of the rational Ej model 

n = (one eigenstate) 

0o = 1 , e = . 

n = 1 (one eigenstate) 

14 

0i = r 2 (I + I81/), ei = 2w. 
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n = 2 (four eigenstate) 

36 252 
<h,i = r| (1 + 14z/)r 2 + — (1 + 14i/)(l + 18u) , e 2 ,i = 4kj , 

UJ UJ 

02.2 = r 6 + -(1 + 10z/)rf - — =-(1 + 10z/)(l + 14z/)r 2 + — =-(1 + 10i/)(l + 14i/)(l + 18i/) , 

£2,2 = 6w , 

c 1 c 1 on 

02.3 = r 8 - -(1 + 6z/)r 6 - -^(1 + 61/) (1 + 10i/)r| + -^-(1 + 6i/)(l + 10*/) (1 + 14i/)ts 

630 

-^(1 + 6i/)(l + 10z/)(l + 14i/)(l + 18i/) , 

£2,3 = 8w , 

02.4 = no + + 2z/)r 8 - (1 + 2i/)(l + 6z/)r 6 - ^3 (1 + 2i/)(l + 61/) (1 + 10i/)r 2 2 

45 

+ 2-i(l + 2i/)(l + 6i/)(l + 10i/)(l + 14!/)r 2 
63 

-^(1 + 2i/)(l + 6i/)(l + 10i/)(l + 14i/)(l + 18*/) , 
e 2 , 4 = lOw , 



First eigenfunctions of the rational E% model 

n = (one eigenstate) 

00 = 1 , eo = 0. 

n = 1 (one eigenstate) 

g 

01 = r 2 - -(l + 30i/) , ei = 2w. 

UJ 

n = 2 (one eigenstate) 



20 _ . x 80 



02 = r 2 2 (1 + 24i/) r 2 + (1 + 24i/)(l + 30i/) , e 2 = 4w. 



47 



n = 3 (two eigenstates) 



q 36. . o 360 , . . . 960 , 
03,1 = r| (l+20i/) T |+— =- (1+20i/)(1+24i/) 75 T (1+20i/)(1+24i/)(1+30i/) , 

<jj 

e 3j i = 6w , 

21 189 
03 ' 2 = T8 + 40^ (1 + 18iy) ^ " 20^ C 1 + 18iy ) C 1 + 20l/ ) r l+ 

^ (1 + 18i/)(l + 20i/)(l + 24*/) r 2 - ^ (1 + 18i/)(l + 20i/)(l + 24i/)(l + 30i/) , 

£3,2 = 8a> . 

n = 4 (two eigenstates) 



ki = r 2 " - ( 7 + 120z/ ) T 2 3 + (1 + 20i/) (7 + 120i/) r|- 



960 1920 

-j3- (1 + 20i/)(l + 24z/)(7 + 120i/) r 2 + (1 + 20i/)(l + 24i/)(l + 30i/)(7 + 120i/) , 



e4 5 i = 8 a; , 

21 24 21 

04,2 = r 2 r 8 + — - (1 + 18z/) r 2 4 (1 + 10i/) r 8 - (1 + 15z/) (1 + 18i/) r|+ 

40 

■^3" (1 + 15i/)(l + 18i/)(l + 20i/) r| - —j- (1 + 15i/)(l + 18i/)(l + 20i/)(l + 24i/) r 2 + 
2016 



a; 



5 



(1 + 15i/)(l + 18i/)(l + 20i/)(l + 24i/)(l + 30z/) , 
£4,2 = 10 a; . 
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